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Executive  Summary 

The  scientific  vision  for  this  EOARD/AFRL  funded  project  is  the  development  of  a 
novel  and  highly  ordered  complex  architecture  fibrous  material  for  additive  layer 
manufacturing.  Using  novel  chemistry  for  the  generation  of  new  hierarchical  fibre 
architectures  and  controlled  field-effect  assembly,  functionally  graded,  stiffness 
modulated  architectures,  analogous  to  those  found  in  nature,  are  synthesised  to 
realise  enhanced  mechanical  performance,  multi-dimensional  composite  structures. 

This  report  provides  the  detailed  experimental  investigations  into  the  generation  of 
Zinc  oxide  (ZnO)  nanorods  grown  in  aqueous  solutions  onto  the  surface  of  E-glass 
short  fibres  (approximately  3  mm  in  length),  and  a  continuous  borosilicate  glass  fibre 
60  mm  in  length.  This  shows  that  the  nanorods  do  not  compromise  the  mechanical 
properties  of  the  fibres  and  furthermore  the  TEOS  coating  could  'heal'  defects  and 
cause  a  slight  increase  in  fibre  cross  sectional  area.  In  the  second  half  of  the  study, 
both  theoretical  analysis  and  experimental  development  were  used  to  evaluate  the 
potential  of  field  gradient  technologies  at  generating  (1)  micron  length  scale  aligned 
reinforcement  and  (2)  nano  length  scale  aligned  composite  materials,  with  regard  to 
rapid  prototyping  manufacture  of  complex  fibre  constructs.  To  achieve  this,  a  new 
type  of  ultrasonic  device  was  developed,  separating  the  acoustic  system  from  the 
resin  cavity  to  allow  for  the  easy  manufacture  of  multiple  samples.  Furthermore,  an 
effective  method  has  been  outlined  to  manufacture  and  mechanically  characterise 
the  micron  length  scale  reinforced  thin  laminar  composites.  With  an  8  %  difference 
in  stiffness  between  the  two  principal  directions,  anisotropy  was  demonstrated  for 
unidirectionally  reinforced  discontinuous  fibre  composites  under  uniaxial  tensile 
loading.  A  43  %  improvement  in  strength  was  observed  for  samples  tested  parallel  to 
the  direction  of 
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the  fibre  reinforcement  over  those  strained  normal  to  the  fibre  direction,  despite  the 
relatively  low  volume  percentage  of  the  reinforcement  phase.  This  technique  shows 
great  potential  for  the  low  cost  instantaneous  alignment  of  structural  reinforcement 
to  generate  the  light-weight  high  performance  structures  required  for  future  military 
requirements. 
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1.  Introduction 

In  nature,  both  material  and  structure  are  formed  according  to  the  principles 
of  biologically  controlled  self-assembly,  a  process  defined  as  the  spontaneous  and 
reversible  ordering  of  small  molecular  building  blocks  under  the  influence  of  non- 
covalent,  static  interactions.  Further,  natural  materials  have  been  shown  to  possess 
outstanding  mechanical  properties  that  are  generally  accredited  to  the  functional 
arrangement  of  short  fibrous  reinforcement  entities  over  many  levels  of  hierarchy. 


Diameter  of 


Figure  1.1:  Simplified  representation  of  the  tendon  structure,  decomposed  into  four  levels 
of  hierarchical  organisation.  Tendon  (a)  is  made  up  of  several  parallel  fasciles  (b),  each 
containing  a  large  number  of  collagen  fibrils  (c),  which  in  turn  represent  assemblies  of 
highly  aligned  staggered  collagen  molecules  (d).  [1] 

There  are  numerous  examples  of  biological  systems  with  varying  degrees  of 
fibrillar  organisation  that  can  be  found  in  nature.  For  the  purposes  of  this  study,  the 
specific  characteristics  of  micron-size  collagen  fibrils  have  been  selected.  Collagen  is 
a  major  constituent  of  tendon  (as  well  as  many  other  biological  tissues)  for  which, 
importantly,  the  following  relation  has  been  proven  to  hold  true. 

£p  £p  £]yj 

In  this  equation,  eT,  £F,  and  £M  signify  the  tensile  strains  of  the  whole  tendon, 
a  single  collagen  fibril,  and  an  individual  collagen  molecule,  respectively.  While  the 
magnitude  of  strain  is  not  independent  of  strain  rate,  it  is  still  clear  that  only 
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hierarchical  effects  can  cause  the  above  relationship  to  hold.  Specifically,  it  has  been 
suggested  that  the  overall  extension  of  the  tendon  may  be  obtained  by  summing 
over  fibril  elongation  and  deformation  in  the  matrix  due  to  shear  [1],  Similarly, 
molecular  strain  and  intermolecular  shear  are  thought  to  constitute  total  fibril  strain 
[1].  Collagen  fibrils  are  assembled  by  an  axial  staggering  of  individual  collagen 
molecules  (Figure  1.1,  d)  that,  when  suspended  in  water,  tend  to  adapt  a  triple  helix 
conformation.  With  each  helix  measuring  approximately  300  nm  in  length  and  1.5 
nm  in  diameter,  molecules  are  arranged  at  a  periodicity  of  67  nm  in  the  longitudinal 
direction,  a  spacing  of  35  nm  between  ends,  and  a  varying  separation  of  1.1  to  1.6 
nm  between  neigbouring  strands  [1,  2],  A  small  number  of  covalent  cross-links  can 
be  detected  between  helices.  Whole  fibrils  are  typically  found  to  have  diameters  of 
the  order  of  several  hundred  nanometres.  Given  the  above  equation,  it  is  postulated 
that  the  architectural  design  features  analogous  to  those  of  collagen  could  offer 
significant  improvements  in  the  overall  mechanical  performance  of  currently 
available  short  fibre  composite  components.  Although  relatively  good  alignment  has 
already  been  achieved  for  randomly  dispersed  particles  in  solution  [3-5],  no  means 
are  presently  accessible  allowing  reinforcement  entities  to  be  ordered  in  accordance 
with  specific  patterns. 

The  goal  of  this  research  is  the  development  of  a  novel  architecture,  high 
performance  fibre  reinforced  polymer  composite  material  incorporating  hierarchical 
mutually  off-set,  but  highly  aligned  fibre  reinforcement  with  integration  and 
directional  control  across  nano,  micro  and  meso  length  scales.  High  performance 
structural  reinforcement  targeted  for  additive  layer  manufacturing  offers  the 
potential  for  'building'  multifunctional  composite  materials  in  both  controlled 
manufacturing  environments  and  at  the  front  line  ensuring  the  timely  introduction 
of  new  designs  in  a  rapidly  changing  theatre  of  conflict. 
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2.  Research  Methodology 

The  work  undertaken  within  this  12-month  Seed  Project  study  is  primarily 
experimental  with  some  supporting  numerical  components.  The  specific  scientific 
objectives  set  out  at  the  onset  of  the  research  project  are  captured  below: 

[1]  To  control  reinforcement  alignment  and  fibre  staggered  spacing  by 
ultrasonic  manipulation. 

[2]  To  select  and  adapt  a  host  polymer  matrix  (e.g.  polyester,  vinylester  or 
epoxy)  to  make  it  suitable  for  suspension  and  'fixing'  of  the  aligned  fibre 
configuration  after  ultrasonic  manipulation. 

[3]  To  maximise  fibre  volume  fraction  with  hierarchical  reinforcement 

[4]  To  experimentally  and  analytically  characterise  the  material  produced  in  [2] 
and  [3], 

To  address  these  objectives  three  parallel  investigations  were  undertaken, 
namely: 

a)  Development  of  hierarchical  fibres  (for  ultrasonic  alignment) 

b)  Ultrasonic  alignment  of  nanowhiskers  reinforcement 

c)  Ultrasonic  alignment  of  micron  length  scale  reinforcement 

The  results  of  these  investigations  are  reported  in  the  following  sections. 
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3.  Ultrasonic  alignment  of  nanowhiskers  reinforcement 

3.1  Introduction 

The  use  of  acoustic  radiation  forces  to  manipulate  particles  has  proven  useful 
across  a  number  of  different  disciplines,  most  notably  within  microfluidics  for 
biomedical  purposes  [1,2]  (e.g  cell  manipulation  and  separation),  using  standing 
waves  with  trap  locations  determined  by  the  dimensions  of  the  system.  There  has 
been  much  less  research  into  the  use  of  such  acoustic  radiation  pressure  to 
manipulate  and  orient  cylindrical  particles  (i.e  fibres),  particularly  from  an 
experimental  perspective.  The  manipulation  and  alignment  of  short  fibres  with  an 
ultrasonic  field  could  provide  a  solution  to  the  problem  of  manufacturing 
unidirectional  short  fibre  composite  materials.  Cellulose  is  an  interesting  material  for 
composite  manufacturing,  as  it  is  abundant  and  cheap,  and  an  environmentally 
friendly  alternative  to  synthetic  fibres.  The  use  of  cellulose  nanowhiskers  could 
produce  composites  with  improved  mechanical  properties  due  to  the  high  tensile 
strength  of  the  whiskers,  which  would  remove  one  of  the  largest  issues  with  current 
cellulose  reinforced  composites.  Due  to  the  small  size  of  the  whiskers,  however, 
conventional  techniques  for  producing  unidirectional  reinforcements  are  ineffective. 
Previous  work  has  shown  the  efficacy  of  using  strong  magnetic  fields  to  align 
cellulose  nanowhiskers[3]  (CNW's)  for  the  purpose  of  producing  nanocomposite  [4], 
It  has  been  shown  that  3-dimensional  and  layered  structures  can  be  formed  using 
magnetic  field  alignment  of  particles  coated  with  superparamagnetic  nanoparticles, 
but  in  order  to  optimise  the  system  for  minimal  magnetic  field  strength  the  size  and 
shape  of  the  particles  are  constrained  [5],  Acoustic  trapping  can  provide  solutions  to 
situations  where  other  trapping  techniques  are  unsuitable,  as  is  commonly  the  case 
with  the  manipulation  of  cells  with  optical  tweezers,  as  the  laser  intensity  required 
for  trapping  can  cause  damage  to  the  cells  [6],  Acoustic  trapping  also  has  the 
advantage  of  being  able  to  manipulate  large  numbers  of  particles  simultaneously 
with  relative  ease,  making  this  an  ideal  technique  for  particle  filtering. 

Much  of  the  previous  work  with  sonotweezers  has  revolved  around  the 
principle  of  manipulating  particles  to  desired  positions  by  varying  the  position  of  the 
trap,  which  has  successfully  been  achieved  by  producing  standing  wave  from  two 
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counter-propagating  waves  and  varying  the  phase  between  the  transducers  [7],  In 
this  work,  the  aim  is  to  produce  a  patterned  distribution  of  nanowhiskers  with  a 
common  orientation,  with  no  need  for  further  manipulation,  thus  allowing  the 
acoustic  field  to  remain  unchanged  throughout  the  process.  The  most  important 
factor  then  becomes  the  use  of  an  appropriate  frequency  to  maximise  the  forces  on 
the  nanowhiskers. 


PMMA 


Figure  3.1:  Schematic  representing  the  transducers  producing  a  2MHz  standing  wave  field 
with  counter  propogating  waves.  The  dashed  lines  show  the  pressure  nodes  where  dense  particles 


are  trapped 


A  schematic  of  the  initial  rig  design  is  shown  in  Figure  3.1.The  transducers  are 
lead  zirconate  titanate  (PZT)  piezoceramic  plates,  and  are  generally  used  at  their 
lowest  frequency  mode  of  2MHz.  The  dimensions  of  the  water  filled  cavity  are  80mm 
x  30mm,  which  is  slightly  larger  than  a  microscope  slide  so  that  samples  can  be 
produced  on  a  slide  and  removed.  Within  Figure  3.1,  the  horizontal  lines  represent 
the  nodal  pressure  positions  where  dense  particles  are  trapped. 

The  frequencies  used  in  this  work  are  unlikely  to  align  individual 
nanowhiskers,  and  therefore  no  surfactant  was  used  to  produce  a  stable  dispersion 
of  separate  whiskers.  Instead,  it  is  hoped  that  by  manipulating  clumps  of 
nanowhiskers  into  traps,  the  concentrations  of  nanowhiskers  within  the  confined 
regions  would  be  large  enough  to  form  a  nematic  phase,  thus  creating  the  desired 
alignment  [8],  The  local  concentrations  required  to  form  such  a  liquid  crystal, 
however,  are  very  high,  and  it  is  not  known  whether  the  acoustic  radiation  forces  on 
the  fibres  will  be  sufficient  for  long  range  rotational  order  to  become  energetically 
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favourable.  It  is  hoped  that  the  fibres  may  partially  align  due  to  the  acoustic  pressure 
gradient  prior  to  confinement  within  the  traps.  There  has  been  some  investigation 
into  the  alignment  and  trapping  processes  separately,  and  a  semi-empirical  model 
has  been  shown  to  describe  the  associated  velocities  of  these  processes  [9],  The 
equaltion  is  given  below: 


?>(:r)  oc  —  rsin[2kx\ 
i o{0)  oc  —  (r/l2)sin[20] 

Here  'v'  is  the  velocity  of  the  centre  of  mass  of  the  fibre  towards  a  trap 
position  and  'co'  is  the  angular  velocity  of  the  fibre  towards  an  aligned  position,  V,  Y 
and  T  are  the  position,  radius  and  length  of  the  fibre  respectively,  and  '0'  is  the 
angle  of  the  fibre  away  from  the  aligned  orientation.  As  the  angular  velocity  of  the 
fibre  is  dependent  upon  the  radius  over  the  length  squared,  the  alignment  process 
will  occur  more  slowly  for  longer  fibres  and  fibres  with  a  large  aspect  ratio. 

Previous  work  has  been  performed  within  ACCIS  to  align  glass  microfibres 
with  acoustic  radiation  pressure  [10].  Within  this  research,  epoxy-glass  composite 
materials  with  ultrasonically  aligned  fibres  have  been  produced  in  thin  panels,  and 
the  resulting  material  exhibit  anisotropic  in-plane  mechanical  properties.  To  date, 
the  aim  has  been  to  produce  unidirectional  composite  materials  with  a  simple 
acoustic  field  produced  by  opposing  flat  transducers.  Theoretically,  acoustic 
radiation  pressure  could  enable  more  complex  3-dimensional  arrangements  of  the 
reinforcing  fibres  to  be  achieved. 

The  initial  theoretical  work  to  determine  the  force  on  particles  due  to 
acoustic  radiation  pressure  is  accredited  to  King  [11],  A  number  of  assumptions  were 
made  in  this  initial  work,  as  the  force  was  determined  for  incompressible  spherical 
particles  in  an  inviscid  medium.  More  recently,  extensions  to  this  work  have  been 
made  which  take  into  account  the  compressibility  of  the  spheres  and  the  viscosity  of 
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the  medium  [12,13],  and  some  which  determine  the  force  on  a  cylindrical  particle 

[14]. 

The  acoustic  pressure  on  spherical  particles  causes  a  force  which  can  be 
considered  as  simply  acting  about  the  centre  of  mass,  resulting  in  a  translational 
displacement  of  the  particle  and  ignoring  any  change  in  angular  momentum.  In  the 
case  of  a  cylindrical  particle,  however,  the  variation  in  the  local  pressure  along  the 
length  of  the  particle  can  result  in  a  rotational  torque  as  well  as  this  translational 
force  [9],  It  is  clear  that  the  dimensions  of  the  particle  determine  the  optimal 
frequency  for  manipulation.  It  has  been  shown  that  the  length  of  the  fibre  must  be 
less  than  a  quarter  of  a  wavelength  [16],  which  is  an  obvious  constraint  to  prevent 
the  fibre  'bridging'  across  neighbouring  pressure  nodes  and  thus  experiencing  no  net 
aligning  force.  However,  in  the  case  that  the  length  of  the  fibre  is  many  orders  of 
magnitude  smaller  than  the  wavelength  of  the  acoustic  field,  the  pressure  difference 
across  the  length  of  the  fibre  will  be  small  due  to  the  low  pressure  gradient.  This  will 
result  in  small  aligning  torques,  which  may  not  be  sufficient  to  even  partially  align 
the  fibres.  It  is  crucial,  therefore,  to  optimise  the  wavelength  of  the  field  for  the 
given  fibre  sizes.  Problems  will  occur  when  the  length  of  the  fibres  is  not  uniform, 
especially  in  the  case  of  unusually  large  fibres,  as  small  fibres  may  not  be  aligned  but 
are  generally  confined  within  a  trap.  Large  fibres  are  not  trapped,  however,  and 
intermolecular  forces  mean  that  smaller  fibres  collect  around  a  large  fibre, 
preventing  their  alignment. 

3.2  Pressure  field  monitoring 

A  model  was  produced  which  mapped  the  predicted  acoustic  field  pressure 
amplitude  between  two  transducers  placed  at  a  given  distance  from  one  another. 
With  the  transducers  perfectly  aligned  and  parallel  to  each  other,  vibrating  at  the 
same  frequency  and  amplitude,  a  simple  standing  wave  pattern  was  produced,  as 
shown  in  Figure.  In  order  to  determine  the  resulting  acoustic  field,  a  finely  spaced 
Cartesian  grid  was  created  such  that  there  were  at  least  10  grid  points  per 
wavelength  to  provide  a  high  resolution  for  plotting  purposes.  This  grid  was  used  as 
the  area  over  which  the  pressure  amplitude  was  calculated  and  finally  plotted.  The 
transducers  were  defined  by  lines  between  two  points,  also  in  Cartesian  coordinates, 
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which  were  then  split  into  individual  Huygens  sources,  ensuring  that  the  spacing 
between  each  point  source  was  less  than  or  equal  to  half  a  wavelength  to  prevent 
the  formation  of  grating  lobes.  The  transducers  were  defined  outside  of  the  grid  area 
to  avoid  singularities  caused  by  zero  distance  values.  The  total  distance  from  each 
transducer  point  source  to  each  point  in  the  grid  was  then  calculated,  and  the 
acoustic  pressure  amplitude  was  calculated  using  the  following  expression: 


Vi  = 


,ikr 


Vr 


Here,  ‘r’  is  the  total  distance  of  the  grid  point  from  the  transducer  source  point, 
‘k’  is  the  wavenumber  and  ‘pi’  is  the  pressure  resulting  from  a  single  Huygens  source 
point.  The  result  of  this  is  essentially  a  circular  wavefield  pattern  from  each  point 
source,  from  which  the  total  pressure  amplitude,  ‘ptotai\  produced  by  the  transducers  is 
computed  through  superposition  by  summing  over  all  of  the  source  points  and  taking 
the  absolute  value  of  this  complex  result: 


Vtotal  — 


I 


Vi 


x  (mm) 


Figure  3.2:  Acoustic  pressure  amplitude  at  the  centre  point  between  two  perfectly  parallel 
transducers  separated  by  80mm.  Both  transducers  are  operating  at  a  frequency  of  2M Hz.  The 
plotted  region  corresponds  to  5  wavelengths  in  the  x-direction  by  13  wavelengths  in  the  y-direction 

According  to  the  theory  discussed  previously  the  use  of  high  frequencies  is 
essential  to  manipulate  individual  CNW's  [9],  as  with  low  frequencies  the  pressure 
gradient  will  not  be  sufficient  to  align  nanoscale  fibres.  As  a  result  of  this,  any 
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misalignment  of  the  transducers  becomes  increasingly  critical,  as  this  causes  a 
greater  phase  change  of  the  wavefronts  along  the  length  of  the  transducer.  The 
effects  of  this  are  shown  in  Figure  3.3.  At  a  frequency  of  2MHz,  the  shape  and 
position  of  the  acoustic  traps  with  both  transducers  rotated  anticlockwise  by  2°  are 
very  similar  to  those  for  perfectly  aligned  transducers  (Figure  3.2).  Flowever,  for  the 
same  2°  at  6.8MFIz,  which  roughly  corresponds  to  the  second  resonant  mode  of  the 
transducers  used,  the  shape  of  the  traps  has  been  significantly  distorted.  The 
maximum  pressure  of  the  system  is  reduced  by  approximately  30%  for  a  2°  rotation 
at  6MFIz,  while  the  same  rotation  at  2MFIz  only  results  in  a  10%  reduction  in  the 
maximum  pressure  of  the  system.  The  minimum  pressure  amplitude  of  the  system  is 
increased  as  a  result  of  the  rotation  of  the  transducers,  with  the  increase  in  the 
6.8MFIz  field  being  3  times  greater  than  in  the  2MFIz  field.  As  the  periodicity  of  the 
wavefields  remains  the  same,  both  of  these  factors  result  in  a  smaller  acoustic 
pressure  gradient,  which  is  more  significantly  affected  with  increasing  frequency. 


x  (mm) 


x(tnm) 


Figure  3.3:  Acoustic  pressure  field  for  the  same  transducers  as  in  Figure  3.2  over  the  same 
spatial  area,  with  both  transducers  rotated  by  2  degrees  anticlockwise.  The  operating  frequencies 

are  6MHz  (top)  and  2MHz  (bottom) 
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3.3  Experimental  Validation 

As  aforementioned,  misalignment  of  the  transducers  can  lead  to  significantly 
reduced  acoustic  radiation  forces.  To  improve  the  transducer  alignment,  a  rig  was 
produced  which  allowed  for  the  transducers  to  be  rotated,  as  shown  in  Figure  3.4. 
This  also  enabled  a  secondary  experiment  to  be  performed,  wherein  the  transducers 
were  deliberately  rotated  through  a  certain  angle,  and  the  locations  of  the  traps  were 
determined  by  tracking  the  position  of  polystyrene  microparticles  in  a  solution  within 
the  rig.  The  results  of  this  were  then  compared  with  the  expected  acoustic  field 
predicted  by  the  model. 


80  mm 


Figure  3.4:  Rig  design  allowing  for  rotation  of  transducers 

To  perform  the  initial  calibration  to  ensure  the  transducers  were  as  ideally  aligned  as 
possible,  one  of  the  transducers  was  connected  to  a  signal  generator,  and  the  other  to  an 
oscilloscope.  The  transducers  were  initially  aligned  by  eye,  and  the  signal  generator  was  then 
used  to  drive  one  transducer  at  its  lowest  resonant  mode,  roughly  2MHz,  at  the  maximum 
possible  voltage,  10VPP  (Volts  peak-to-peak).  Higher  voltages  are  possible  with  the  use  of  an 
amplifier  but  generally  lead  to  acoustic  streaming.  By  rotating  the  transducers  slowly,  the 
oscilloscope  was  used  to  determine  the  point  at  which  the  signal  response  had  maximum 
amplitude,  which  corresponds  with  parallel  alignment  of  the  transducers 

In  order  to  verify  the  model,  a  solution  of  1pm  diameter  polystyrene  spheres  was 
placed  in  the  rig.  These  particles  are  easily  trapped  by  the  acoustic  radiation  pressure,  and 
with  suitable  concentrations  they  provide  a  clear  indication  of  the  locations  of  the  traps  with 
a  given  acoustic  field.  200pl  of  10%  polystyrene  microsphere  solution  was  mixed  with  4ml  of 
water  and  placed  within  the  rig  cavity.  The  transducers  were  calibrated  to  be  aligned  and 
driven  with  30VPP  at  a  frequency  2MHz  until  the  trap  sites  were  clearly  visible  by  the 
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locations  of  the  microspheres,  which  took  roughly  1  minute.  Following,  the  transducers  were 
rotated  by  2°  anticlockwise  and  then  driven  at  the  same  frequency  and  voltage  to  compare 
the  effects  of  misalignment  with  those  predicted  by  the  model.  The  results  of  this  are  shown 
in  Figure  3.5. 


Figure  3.5:  lpm  polystyrene  microsplieres  trapped  in  a  2MHz  acoustic  field  (a)  with  aligned 
transducers  (b)  with  transducers  rotated  by  2° 

It  can  be  seen  in  Figure  3.5  that  with  aligned  transducers,  the  trapping  areas 
are  reasonably  uniform  straight  lines,  with  some  waviness.  With  a  2°  rotation  of  the 
transducers  however,  the  traps  are  also  slightly  rotated,  and  the  edges  of  the  traps  are 
significantly  distorted,  as  expected  in  Figure  4.  Overall,  however,  the  traps  remain 
distinct  and  the  shape  is  not  globally  affected.  The  same  procedure  was  performed 
with  a  frequency  of  6.77MHz  and  a  driving  voltage  of  4Vpp,  and  the  results  are 
shown  in  Figure  3.6. 
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Figure  3.6:  lpm  polystyrene  microspheres  trapped  in  a  6MHz  acoustic  field  (a)  with  aligned 
transducers  (b)  with  transducers  rotated  by  2° 

The  driving  voltage  used  at  6.77MHz  was  significantly  lower  than  with  2MHz 
to  minimise  acoustic  streaming  whilst  still  applying  sufficient  forces  to  the 
microspheres  for  manipulation.  As  a  result  the  particles  took  roughly  5  minutes  to 
move  to  the  pressure  nodes.  With  aligned  transducers,  the  traps  form  straight  and 
parallel  lines,  as  predicted  by  the  model.  With  transducers  rotated  by  2°,  the 
effective  trapping  area  was  dramatically  reduced,  with  only  a  very  narrow  band 
forming  an  effective  trapping  site.  This  agrees  well  with  what  is  predicted  by  the 
model. 
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3.4  Experimental  Validation 

3.4.1  Cellulose  Nano  whisker  Alignment 

The  initial  calibration  procedure  was  performed  to  align  the  transducers.  To 
determine  whether  the  traps  were  clearly  defined  parallel  lines  as  expected,  a  glass 
microfibre  solution  was  placed  within  the  rig,  and  their  alignment  within  the 
pressure  traps  was  checked.  Once  this  was  completed,  the  solution  was  removed 
and  the  rig  cleaned  in  preparation  for  using  a  CNW  solution. 

The  cellulose  nanowhiskers  used  were  provided  by  Professor  Stephen 
Eichhorn  from  the  CEMPS  Department  at  Exeter  University.  The  nanowhiskers 
supplied  are  in  a  dry  form,  and  in  this  state  they  are  in  large  clumps.  Two  different 
samples  were  made,  wherein  these  clumps  were  broken  down  to  varying  degrees. 
Whilst  the  specific  batch  of  nanowhiskers  used  in  this  project  have  not  been 
characterised,  similar  nanowhiskers  produced  with  the  same  method  by  the  same 
group  have  been  [8],  These  nanowhiskers  are  derived  from  cotton,  and  have 
dimensions  of  around  130nm  x  20nm,  giving  them  a  relatively  small  aspect  ratio  of 
around  6.  To  determine  whether  individual  nanowhiskers  (or  at  least  very  small 
clumps)  could  be  manipulated,  a  solution  of  the  nanowhiskers  was  produced  by 
mixing  the  nanowhiskers  in  water  and  then  sonicating  for  60  minutes.  The  second 
sample  contained  larger  clumps  of  the  nanowhiskers,  and  was  produced  by  mixing 
the  nanowhiskers  in  water  and  stirring  to  break  down  the  larger  clumps.  In  both 
cases,  only  fibres  from  near  the  surface  of  the  solution  were  used  to  produce  the 
sample  to  avoid  using  excessively  large  clumps  of  the  nanowhiskers  which  settled  at 
the  bottom  of  the  solution. 

A  microscope  slide  was  placed  on  the  floor  of  the  rig  cavity,  so  that  dried 
samples  of  aligned  nanowhiskers  could  be  removed  and  the  rig  reused.  The  rig  was 
filled  with  the  nanowhiskers  solution  such  that  the  transducers  were  just  fully 
submerged.  The  signal  generator  was  then  used  to  apply  a  standing  wave  field  with  a 
frequency  of  2MHz  (giving  a  wavelength  of  736pm)  using  a  driving  voltage  of  lOVpp 
.At  this  stage  the  alignment  process  occurs,  with  the  nanowhiskers  sitting  at  the 
pressure  nodes  of  the  acoustic  field.  A  needle  syringe  is  then  used  to  remove  water 
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slowly  from  the  corner  of  the  rig  to  reduce  drying  time.  The  acoustic  pressure  field  is 
applied  for  30  minutes,  after  which  point  the  signal  generator  is  switched  off.  The 
sample  is  then  left  undisturbed  to  dry  out  over  a  period  of  24  hours,  at  which  point  a 
dry  sample  of  lines  of  nanowhiskers  is  left  on  the  microscope  slide  in  the  bottom  of 
the  rig,  which  can  then  be  removed  and  analysed. 

To  determine  whether  any  alignment  of  the  nanowhiskers  was  present. 
Polarised  Raman  Spectroscopy  was  used.  Monochromatic,  polarised  laser  light  is 
focused  on  a  small  area  of  the  sample  (2  micron  diameter).  This  light  then  excites 
vibrational  modes  within  the  system  into  virtual  energy  states.  When  the  vibrational 
modes  then  relax  to  a  lower  energy  mode,  they  emit  a  photon  with  the  same  energy 
as  the  energy  difference  between  the  initial  excited  mode  and  the  relaxed  state. 
Certain  vibrational  modes  are  polar,  such  that  the  polarisation  of  the  incident  laser 
light  must  align  with  the  relevant  vibration  of  the  molecular  structure  in  order  to 
cause  an  excitation.  Such  an  example  exists  within  the  cellulose  nanowhiskers  used 
in  this  case.  Along  the  cellulose  backbone,  vibration  of  the  C-O-C  bond  results  in  a 
peak  at  a  Raman  shift  of  1095cm-l,  and  the  intensity  of  this  peak  is  greatly  increased 
when  the  polarisation  of  the  incident  light  is  along  the  length  of  the  cellulose  chain. 
In  the  case  of  cellulose  nanowhiskers,  this  corresponds  to  the  long  axis  of  the 
whisker.  Therefore,  if  there  is  overall  alignment  of  the  nanowhiskers,  the  intensity  of 
this  peak  should  be  highest  when  along  the  direction  of  alignment.  In  order  to 
determine  whether  there  was  any  alignment,  Raman  spectra  were  taken  of  a  point 
along  one  of  the  'stripes'  of  nanowhiskers  of  a  dried  sample  while  the  laser 
polarisation  was  rotated,  with  one  spectrum  being  taken  every  3°  from  0°  to  90°.  The 
spectra  were  then  fitted  to  one  another  by  removing  signal  noise  (as  a  result  of 
fluorescence)  and  correcting  for  variations  in  the  incident  laser  intensity,  in  order  to 
determine  the  real  variation  in  the  intensity  of  the  peak  at  1095cm-l  and  therefore 
the  degree  of  alignment  of  the  sample.  A  second  characteristic  peak  of  cellulose 
occurs  at  a  Raman  shift  of  around  2900cm-l ,  which  results  from  the  vibration  of  the 
C-H  bonds  within  cellulose.  These  vibrations  are  non-polar,  such  that  the  polarisation 
of  incident  laser  light  has  a  negligible  effect  on  the  intensity  of  the  Raman  peak.  By 
taking  the  ratio  of  the  intensity  of  the  peaks  at  1095cm-l  and  2900cm-l ,  this  gives  a 
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value  proportional  to  the  fraction  of  interacting  cellulose  that  is  aligned  with  the 
polarisation  of  the  laser.  This  significantly  simplifies  the  spectral  fitting  procedure, 
and  reduces  errors  resulting  from  any  variation  in  laser  intensity  due  to  rotation  of 
the  polariser. 

3.4.1  Frequency  Modulation 

One  of  the  primary  issues  with  aligning  fibres  within  an  acoustic  pressure 
field  occurs  as  a  result  of  the  short  timescale  over  which  the  alignment  process 
occurs,  particularly  in  a  low  viscosity  medium  such  as  water.  As  a  result  of  this,  the 
fibres  move  swiftly  into  the  traps,  with  the  frictional  forces  between  fibres  tending 
to  resist  the  alignment  process,  resulting  in  randomly  aligned  clumps  of 
nanowhiskers  confined  within  the  traps.  In  order  to  counter  this,  a  method  was 
devised  to  repeatedly  'relax'  the  fibres  away  from  the  trap  and  then  force  them  back 
to  the  traps.  To  achieve  this,  the  frequency  of  the  transducers  was  modulated 
between  2MHz  and  6MHz.  The  lower  frequency  mode  produced  a  larger  force  on 
the  fibres,  and  so  they  were  gradually  confined  within  the  trap  sites  produced  by  this 
frequency. 

3.5  Experimental  Results  and  Discussion 

3.5.1  CNW  Alignment 

The  results  of  the  initial  calibration  using  glass  fibres  shows  that  the  traps  are 
well  defined  and  align  fibres  on  this  length  scale  effectively. 

Two  samples  were  analysed  using  the  Raman  spectroscopy  technique 
described  previously.  Both  samples  were  produced  with  a  2MHz  applied  acoustic 
field,  and  allowed  to  dry  at  room  temperature.  An  example  of  the  CNW  solution  with 
the  acoustic  field  applied  is  shown  in  Figure  3.7. 


R.  S.  Trask  @  Bristol  ac.  uk 


19 


30/03/2014 


Distribution  A:  Approved  for  public  release;  distribution  is  unlimited. 


Additive  Layer  Manufacturing  Of  Biologically  Inspired  Short  Fibre  Reinforced  Composites 

Prepared  by  Dr  R.S.  Trask 


Figure  7:  Optical  Microscopy  (Zeiss  Axio  Imager  M2,  5x  objective)  image  of  CNW's  with  2MHz 

standing  wave  field  applied 

The  spectra  need  correction  for  varying  levels  of  noise  (resulting  from 
fluorescence)  before  the  degree  of  alignment  of  the  fibres  can  be  determined.  A 
baseline  correction  for  the  spectra  was  performed.  It  was  then  assumed  that  the 
intensity  of  the  peak  at  2900cm"1  was  invariant,  which  is  a  reasonable  assumption  if 
the  region  of  the  sample  analysed  does  not  vary  [17].  Having  scaled  the  spectra 
accordingly  to  normalise  against  this  peak,  the  variation  in  the  intensity  of  the  peak 
at  1095cm"1  could  be  determined.  The  corrected  and  normalised  spectra  for  each 
sample  are  shown  in  Figure  3.8.  The  peaks  visible  in  sample  1  at  around  3100cm"1  are 
the  result  of  high  energy  cosmic  rays  incident  on  the  CCD  array. 
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Corrected  Raman  Spectra:  Sample  1 


Corrected  Raman  Spectra:  Sample  2 


Figure  3.8:  Corrected,  normalised  Raman  spectra  of  the  two  dried  samples  of  aligned  CNW's 

To  clearly  visualise  the  variation  of  the  intensity  of  the  peak  with  angle,  a 
polar  plot  is  used,  with  the  angle  between  the  light  polarisation  and  the  axis  of  the 
acoustic  traps  represented  by  0  and  the  intensity  of  the  peak  for  that  spectra 
represented  by  V.  The  results  are  shown  in  Figure  3.9.  As  the  distance  of  each  point 
from  the  origin  represents  the  ratio  between  the  intensity  of  two  Raman  shift  peaks, 
the  value  is  dimensionless,  and  only  the  relative  value  of  different  points  at  different 
angles  can  be  used.  This  gives  a  qualitative  idea  of  the  degree  of  alignment  of  the 
nanowhiskers  in  the  sample,  with  a  perfectly  circular  plot  showing  entirely  random 
orientation  of  the  local  particles  and  a  thin  flat  line  representing  perfect  alignment  in 
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a  single  direction.  As  spectra  were  only  taken  from  0°  to  90°,  these  values  were  then 
mirrored  to  form  a  full  polar  plot  which  gives  a  clearer  impression  of  the  anisotropy 
of  the  sample. 


Polarisation  Angle  vs  Peak  Intensity:  Sample  1 

90 

0.6 


Polarisation  Angle  vs  Peak  Intensity:  Sample  2 
60 

1.5 


Figure  3.9:  Variation  in  1095cm1  peak  intensity  with  polarisation  angle. 

The  polar  plots  for  the  two  samples  give  contrasting  ideas  about  the  possible 
alignment  of  the  nanowhiskers  within  the  sample.  Sample  2  appears  to  show  some 
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anisotropy,  with  a  preferred  orientation  along  the  length  of  the  acoustic  traps, 
although  the  maximum  variation  in  intensity  with  angle  is  small  (approximately  4/3). 
The  results  of  sample  1  appear  to  show  a  more  uniform  intensity  with  angle,  and 
therefore  a  lack  of  any  alignment  within  the  sample.  The  spectra  from  sample  1 
generally  had  higher  levels  of  background  noise,  making  analysis  of  the  spectral 
peaks  more  difficult  and  more  prone  to  error,  but  overall  there  is  no  strong 
indication  of  alignment  of  individual  nanowhiskers  within  either  sample.  The  spectra 
for  sample  1  with  polarisation  angles  from  81°  to  90°  had  exceptionally  high 
fluorescence  levels,  which  made  it  impossible  to  accurately  determine  the  intensity 
of  any  peaks  present,  and  so  this  data  has  been  omitted.  Possible  explanations  for 
the  anisotropy  observed  in  sample  2  are: 

•  Acoustic  radiation  pressure  creating  high  concentration  regions  of  CNW's, 
which  self-assemble  into  a  nematic  phase. 

•  Acoustic  manipulation  of  locally  aligned  clumps  of  CNW's  into  trap  sites. 

•  Coincidental  local  alignment  of  sampled  region. 

It  is  unlikely  that  the  observed  anisotropy  is  a  result  of  coincidental  local 
alignment  of  the  whiskers,  as  the  sampling  region  will  have  contained  thousands  of 
whiskers.  Due  to  the  relatively  low  forces  created  by  the  acoustic  radiation  pressure, 
it  is  suggested  that  the  observed  anisotropy  is  a  result  of  the  manipulation  of 
internally  self-oriented  clumps  of  CNW's.  A  representation  of  this  distribution  is 
shown  in  Figure  3.10. 
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Figure  10:  Possible  arrangement  of  clumps  of  aligned  nano  whiskers  which  are  manipulated  as  a 

group  by  the  acoustic  radiation  pressure. 


3.6  Frequency  Modulation 

The  effect  of  switching  between  two  separate  resonant  modes  was  made 
extremely  clear  when  using  the  glass  microfibres.  Using  a  low  modulation  frequency, 
it  was  possible  to  observe  the  movement  of  the  particles  in  real  time  between  the 
adjacent  traps.  A  secondary  effect  of  this  method  meant  that  traps  with  a  high  fibre 
concentration  'leaked'  fibres  to  neighbouring  traps,  producing  a  more  uniform 
distribution  of  fibres  in  each  band,  as  shown  in  Figure  3.11. 


Figure  3.11:  Optical  Microscopy  image  of  glass  microfibre  solution;  (a)  Randomly  oriented  with  no 
acoustic  field  applied,  (b)  aligned  with  single  2M  Hz  frequency,  (c)  aligned  with  frequency 
modulating  between  2MHz  and  6MHz  at  a  frequency  of  10Hz 


There  was  less  obvious  improvement  in  the  case  of  the  CNW's  when  using 
this  frequency  modulation,  which  is  potentially  a  result  of  the  nanowhiskers 
propensity  to  agglomerate  into  clumps  in  solution,  resulting  in  intermolecular  forces 
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which  are  greater  than  the  separating  forces  applied  by  the  modulating  radiation 
pressure.  Despite  this,  it  is  clear  that  for  larger  fibres  with  weaker  intermolecular 
forces  this  modification  of  the  alignment  process  gives  improved  alignment,  albeit  at 
a  slower  rate  than  with  a  single  frequency,  and  shall  be  investigated  in  greater  detail 
during  future  work. 

3.7  Conclusions 

The  results  indicate  that  the  use  of  acoustic  radiation  pressure  is  not  effective 
in  driving  the  formation  of  a  fully  nematic  phase  in  CNW's,  with  only  a  small  degree 
of  orientational  order.  The  rapid  manipulation  of  particles  into  trap  sites  does  not 
seem  to  allow  for  the  whiskers  to  align  initially,  and  the  acoustic  forces  are  then  not 
sufficiently  strong  to  produce  the  concentrations  required  to  form  a  liquid  crystal.  It 
has  been  shown  that  agglomerates  of  nanowhiskers  can  be  manipulated  into  a 
macroscopic  structure  with  partial  alignment,  however,  and  with  more  precise 
equipment  it  is  entirely  possible  that  global  orientational  order  of  individual 
nanowhiskers  may  be  achieved.  As  previously  discussed,  this  would  require  the  use 
of  significantly  higher  frequencies  than  those  implemented  in  this  work,  and  so  the 
effects  of  streaming  may  become  increasingly  important.  The  use  of  simple 
frequency  modulation  has  improved  the  alignment  possible,  indicating  that  there  is  a 
still  much  work  to  be  done  to  optimise  the  process,  and  potentially  reduce  the 
precision  required  in  the  apparatus  to  work  at  higher  frequencies. 

3.8  Future  Work 

Having  demonstrated  the  ability  to  produce  macroscopic  structures  of  CNW's 
within  a  fluid,  it  would  be  desirable  to  produce  a  sample  of  these  CNW's  in  a  polymer 
matrix  for  tensile  testing,  in  order  to  determine  whether  the  reinforcement  is 
anisotropic.  While  the  Raman  results  indicated  that  there  is  only  partial  alignment  of 
the  nanowhiskers,  the  analysis  was  only  performed  on  small  areas  of  the  samples. 
Though  the  samples  would  not  be  expected  to  exhibit  highly  anisotropic  material 
properties,  it  is  important  to  definitively  prove  this  experimentally. 


R.  S.  Trask  @  Bristol  ac.  uk 


25 


30/03/2014 


Distribution  A:  Approved  for  public  release;  distribution  is  unlimited. 


Additive  Layer  Manufacturing  Of  Biologically  Inspired  Short  Fibre  Reinforced  Composites 

Prepared  by  Dr  R.S.  Trask 


The  use  of  ultrasonic  fields  has  proven  to  be  an  effective  method  for  aligning 
fibres,  assuming  that  the  frequency  used  is  appropriate  for  the  fibre  dimensions.  The 
use  of  nanofibres  introduces  a  large  number  of  new  issues,  which  are  largely  down 
to  the  imprecision  of  the  rig  production.  In  order  to  overcome  these  issues,  it  is 
practical  to  make  use  of  existing  manufacturing  techniques  which  are  inherently 
extremely  precise.  Lithographic  techniques  with  nanoscale  precision  are  routinely 
employed  in  the  production  of  Surface  Acoustic  Wave  (SAW)  devices,  by  etching  and 
then  printing  microscale  or  even  nanoscale  interdigital  transducers  onto  a 
piezoelectric  substrate  (commonly  Lithium  Niobate).  The  spacing  between  two 
opposing  digits  corresponds  to  X/2  of  the  surface  wave  produced  (the  spacing 
between  two  digits  on  the  same  side  is  equal  to  X),  as  represented  in  Figure  12  [18]. 
As  a  result  of  using  lithographic  techniques,  SAW  devices  can  be  produced  with 
extremely  high  Q.  values  (>80000)  which  vibrate  at  extremely  high  frequencies 
(>lGHz)  [19],  It  is  a  future  goal  to  repeat  the  process  for  alignment  presented  in  this 
report  using  a  SAW  device  in  order  to  determine  the  alignment  possible  with  more 
precisely  manufactured  equipment.  This  would  effectively  be  an  optimisation  stage, 
with  the  goal  of  creating  an  effective  manufacturing  process. 

There  were  a  number  of  reasons  for  using  cellulose  nanowhiskers  as  the 
reinforcement  in  this  work.  The  use  of  short  fibres  in  composite  materials  mirrors 
materials  seen  in  nature,  and  short  cellulose  fibres  are  present  in  the  wall  of  plant 
cells.  This  structure  contains  concentric  ring  layers,  each  of  which  have  a  different 
fibre  orientation,  much  like  the  composite  laminates  used  commercially.  Little  is 
known  about  the  properties  of  these  nanowhiskers,  however,  making  it  difficult  to 
optimise  the  system  (e.g.  fluid  medium  density/viscosity,  substrate  material)  to 
maximise  the  effect  of  the  acoustic  radiation  pressure.  The  aim  in  future  would  be  to 
use  other  nanomaterials  which  have  more  suitable  acoustic  properties.  Carbon 
nanotubes  are  theoretically  an  excellent  material  for  this  process  due  to  their 
density  and  rigidity,  but  it  is  generally  difficult  to  produce  a  stable  dispersion  of 
individual  nanotubes  which  would  be  required 
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Figure  12:  Overview  of  Surface  Acoustic  Wave  device 

In  future  work,  there  is  to  be  a  greater  focus  on  the  use  of  microfibres  and 
also  macrofibres,  with  the  goal  of  extending  this  process  into  a  continuous 
manufacturing  process,  either  within  existing  Additive  Layer  methods  or  as  an 
entirely  independent  process.  The  promising  results  from  using  magnetic  fields  to 
align  particles  in  a  plane  [4,5]  indicate  that  using  magnetic  and  acoustic  fields  in 
conjunction  could  lead  to  the  ability  to  produce  hierarchical  structures  with  control 
of  reinforcing  particles  on  multiple  length  scales.  Magnetic  fields  can  be  used  to 
produce  a  unidirectional  system  of  fibres,  and  the  acoustic  field  used  to  produce  a 
structure  to  the  position  of  the  fibres  within  the  magnetically  trapped  plane.  This 
may  allow  for  nanofibres  to  be  aligned,  but  optimisation  of  the  concentration  of 
paramagnetic  nanoparticles  would  be  required  to  allow  for  alignment  with 
significantly  altering  the  geometry  and  mechanical  properties  of  the  nanofibres 
themselves. 
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4.  Development  of  hierarchical  fibres  (for  ultrasonic 
alignment). 

4.1  Introduction 

Structural  hierarchy  at  the  macro  scale  is  generally  well  established,  providing 
a  structure  that  constitutes  several  length  scales  such  as  the  Eiffel  tower  [1], 
However,  hierarchy  has  had  limited  use  in  man-made  materials  [2];  engineering 
materials  such  as  metals  and  composites  have  hierarchy  at  a  length  scale  much 
larger  than  observed  in  nature.  Hierarchy  in  a  material  is  sometimes  conjectured  to 
improve  specific  strength  and  stiffness  [1],  allowing  relatively  weak,  more  widely 
available  constituents  to  be  used  to  obtain  high-performance  macro-scale  properties 
[3].  Biological  materials  such  as  bone  are  naturally  occurring  composites  that  exhibit 
a  rich  hierarchical  structure  [1,  4,  5].  Approximately  seven  levels  of  hierarchy  are 
present:  bone  tissue  typically  of  a  compact  or  spongy  nature  (50  cm),  osteons  and 
haversian  canals  (100  mm),  fibre  patterning  and  alignment  (50  mm),  fibril  arrays 
constituting  the  fibre  (10  mm),  mineralisation  of  collagen  fibrils  (1  mm),  formed 
tropocollagen  proteins  (300  nm)  which  are  made  up  of  amino  acids  (1  nm)  [5],  Each 
level  of  hierarchy  has  been  associated  with  various  physical  properties;  for  instance, 
fibril  crack  bridging  at  the  micro-scale  is  one  mechanism  for  increasing  toughness, 
and  stiff  mineral  platelets  at  the  nano-scale  deform  in  tension  to  transfer  stress 
between  discontinuous  mineral  platelets  providing  gracious  failure.  By  using 
principles  of  such  biological  materials  with  current  accessible  synthetic  systems,  this 
may  enable  the  creation  of  the  next  generation  of  man-made  composites  with 
improved  properties  and  functionality  [6], 

Unlike  natural  composites,  current  man-made  composites  have  a  reinforcing 
fibre  that  is  typically  continuous  throughout  a  component  in  order  to  have 
uninterrupted  load  transfer  [7].  However,  this  can  often  result  in  limited 
manufacturability  [8,  9]  and  a  loss  in  the  number  of  levels  of  hierarchy.  As  such, 
strength  and  ductility  in  current  composite  materials,  such  as  glass  or  carbon  fibre 
epoxy,  are  usually  mutually  exclusive  [10].  Current  composites  therefore  tend  to  lack 
gracious  failure  and  perform  in  a  brittle  manner  leading  to  sudden  in-service  failure. 
In  engineering,  this  often  means  overdesign  of  components  and  an  increased  safety 
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factor,  increasing  mass  and  cost  of  the  components.  Discontinuous  man-made 
composites  on  the  other  hand  have  been  shown  to  increase  ductility,  particularly  for 
brittle  matrices  such  as  epoxy  or  polyester  [9],  However,  in-plane  properties 
matching  that  of  continuous  fibre  composites  have  yet  to  be  achieved. 


Figure  4.2:  Carbon  nanotubes  grafted  onto  a  carbon  fibre  surface  through  a  CVD 
process,  a)  Bare  fibre,  b)  CNT  coated  fibre  [12]. 

Considerable  efforts  have  been  made  to  increase  the  level  of  hierarchy  in 
composites  through  the  use  of  nanomaterials  alongside  typical  micro-scale 
reinforcements,  with  particular  interests  in  carbon  nanotubes  (CNT's)  due  to  their 
intrinsic  high  strength  and  stiffness  [11],  When  CNT's  have  been  incorporated  into  a 
matrix  several  properties  have  been  shown  to  improve.  However,  clear  improvement 
of  in-plane  mechanical  properties  has  yet  to  be  demonstrated  [3],  Current 
mechanical  limitations  of  CNT's  include  the  relatively  poor  quality  of  supply, 
difficulties  in  their  dispersion  and  alignment  within  the  matrix,  and  poor  bonding 
between  the  CNT's  and  the  matrix.  CNT's  however  have  been  shown  to  provide 
other  functional  advantages  such  as  electrical  conductivity,  thermal  conductivity, 
increasing  service  temperature  and  flame  retardancy  [3],  To  solve  some  of  the 
mechanical  limitations,  interest  has  grown  in  incorporating  CNT's  via  grafting  to  the 
fibre  surface  through  a  chemical  vapour  deposition  (CVD)  process,  before 
embedding  them  in  a  polymer  [12]  as  shown  in  Figure  4.1.  By  grafting  the  nano¬ 
scale  reinforcement  onto  the  surface  of  a  fibre,  alleviation  of  agglomeration  and 
alignment  issues  becomes  possible,  whilst  improvements  in  through  thickness 
properties  of  a  composite  have  also  been  demonstrated  [13],  It  is  well  known  that 
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the  performance  of  a  composite  not  only  relies  on  the  properties  of  the  individual 
constituents,  but  also  the  interfacial  properties  between  the  fibre  and  the  matrix 

[14] .  Alteration  of  the  surface  of  the  fibre  can  increase  load  transfer  between  the 
matrix  and  the  fibre,  either  through  bonding  or  mechanical  interlocking;  however, 
this  is  not  always  desirable.  Increased  load  transfer  may  induce  larger  amounts  of 
fibre  fracture  in  a  brittle  fibre  /  brittle  matrix  combination,  and  hence  yield 
brittleness  in  the  material,  a  property  that  well  utilised  carbon  fibre  /  epoxy  resins 
already  exhibit.  Hitherto,  a  relatively  ductile  fibre  coating  with  enhanced  mechanical 
interlocking  remains  unknown. 

Zinc  Oxide  (ZnO)  nanorods  are  a  promising  candidate  due  to  their  intrinsic 
ductility  compared  with  CNT's,  with  a  measured  stiffness  of  approximately  29±8  GPa 

[15] ,  which  is  still  greater  than  that  of  typical  epoxy  of  3-4  GPa.  The  exact  value  of 
the  modulus  however  is  still  disputed  and  has  been  shown  to  vary  depending  upon 
the  nanorod  dimensions  [16,  17].  The  relative  ease  and  variety  of  grafting  /  growth 
methods  is  also  a  benefit  for  production  purposes.  It  is  common  to  have  a  two  stage 
process,  first  seeding  the  surface,  followed  by  growth  of  rods  [18].  Growth  can  be 
performed  hydrothermally  in  a  non-hazardous,  low  temperature  environment.  There 
has  been  profound  interest  in  ZnO  nanorods  of  late  mainly  due  to  a  large  range  of 
high  technology  electronic  applications  such  as  field  effect  transistors  [19],  flexible 
electronics  [20],  light  emitting  diodes  [21],  energy  scavenging  [22],  humidity  sensing 
[23],  and  solar  cells  [24],  These  would  offer  fertile  ground  for  exploiting  added 
functionalities  to  composites  by  incorporating  ZnO  nanorods.  Several  morphologies 
of  nanorods  have  been  documented  [25-28],  showing  a  plethora  of  rod  density  and 
alignment.  We  thus  postulate  that  ZnO  nanorod-coated  fibres,  when  incorporated  in 
a  composite,  facilitate  mechanisms  for  not  only  graceful  mechanical  failure,  but  also 
improved  load  transfer  between  discontinuous  fibre  composites.  In  this  work,  ZnO 
nanorods  have  therefore  been  hydrothermally  grown  onto  flat  glass  surfaces  in 
order  to  optimise  the  method  for  controlling  the  rod  morphology.  Subsequently, 
nanorod  growth  on  a  continuous  (long)  and  then  a  discontinuous  (short)  glass  fibre 
surface  has  also  been  performed.  We  then  investigate,  using  an  atomic  force 
microscopy  (AFM)  and  friction  measurement  capability  at  Bristol  Chemistry, 
interfacial  shear  properties  between  nanorod-coated  surfaces.  These  results  provide 
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the  first  concrete  data  for  us  to  critically  assess  the  prospect  of  nanorod-coated 
fibres  for  improved  ductility  in  composite  materials. 

4.1.2  Discontinuous  fibre  composites 

Discontinuous  composites  are  rarely  used  for  primary  structure  components, 
but  more  so  for  secondary  structures  due  to  their  ease  of  manufacture  [8,  9],  The 
processing  methods  often  reduce  manufacturing  time  scales  whilst  also  providing 
potential  for  novel  manufacturing  techniques  such  as  additive  layer  manufacturing 
(ALM)  [29],  Discontinuous  composites  do  however  tend  to  suffer  from  inferior 
stiffness  and  tensile  strength,  often  caused  through  poor  fibre  alignment,  fibre 
lengths  below  the  critical  fibre  length,  Ic,  and  the  difficulty  in  obtaining  high  fibre 
volume  fractions  [30],  The  most  widely  used  model  to  describe  stress  transfer  in  a 
short  fibre  composite  is  the  shear  lag  model,  focusing  on  tensile  stress  transfer 
between  the  fibre  and  matrix  through  interfacial  shear  stress  (IFSS)  [14].  This  leads 
to  the  critical  fibre  length,  Ic,  which  is  the  minimum  length  of  fibre  required  for  the 
fibre  to  fracture  through  IFSS  transfer  [14]. 

In  general  a  high  fibre  content  is  required  in  order  to  achieve  a  higher 
performing  short  fibre  composite  [31],  and  increasing  volume  fraction  (Vf)  has  been 
shown  to  linearly  increase  fracture  toughness  [9].  Increasing  fibre  content  has 
however  been  found  to  reduce  failure  strain,  with  failure  initiation  often  occurring  at 
fibre  ends  and  propagating  across  the  fibre/matrix  interface  and  finally  through  the 
matrix.  Modelling  has  revealed  that  fibre  ends  act  as  stress  concentrations  with  a 
factor  of  10  or  higher  [31],  relating  failure  to  the  number  of  fibre  ends,  with  stress 
often  relieved  through  matrix  flow,  interfacial  debonding,  or  matrix  failure.  The 
interface  is  dependent  not  only  on  the  chemical  bonding  but  also  the  surface 
topography,  and  should  be  controlled  in  a  manner  very  specific  to  a  particular 
fibre/matrix  combination.  It  is  frequently  the  determining  factor  for  the  type  of 
fracture  mechanism  observed  [9],  By  altering  the  interface  with  the  application  of 
ZnO  nanorods,  enhanced  load  transfer  between  discontinuous  fibres  may  be 
possible,  whilst  the  alteration  of  crack  propagation  and  increased  fracture  toughness 
also  has  potential. 
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4.1.3  Grafting  of  CNT's  to  a  fibre  surface 

Grafting  CNT's  onto  the  fibre  surface  has  been  shown  to  improve  mechanical 
properties  such  as  IFSS  by  approximately  15%  [12];  however,  fibre  tensile  strength  is 
often  reduced  by  approximately  15-55%  [3,  32],  This  is  primarily  due  to  the  chemical 
vapour  deposition  (CVD)  method  used  for  growth,  where  a  metallic  catalyst  is  used 
to  'seed'  the  surface  for  CNT  growth.  It  is  thought  that  this  seeding  method  damages 
the  surface,  effectively  etching  away  the  fibre  surface  resulting  in  surface  pitting 
which  causes  a  micro-scale  stress  concentration,  leading  to  lower  tensile  strength 
[3],  Wettability  of  the  CNT  forests  on  the  surface  has  also  been  a  concern,  with 
morphology  of  the  CNT  forest  and  the  viscosity  of  the  resin  having  a  large  impact.  It 
has  been  shown  that  CNT's  with  a  length  of  30  pm  grown  in  square  arrays  had 
difficulty  in  wetting  when  using  an  epoxy  resin  system.  This  caused  large  contraction 
of  the  square  arrays  and  large  central  voids,  however,  when  a  lower  viscosity  was 
used  with  a  similar  CNT  array  no  signs  of  voidage  were  present  [33],  A  measurement 
of  contact  angle  has  also  been  made  using  a  lower  viscosity  resin  and  a  CNT  length  of 
approximately  2-10  pm.  This  showed  a  reduction  in  contact  angle  with  poly(methyl 
methacrylate)  (PMMA)  polymer,  indicating  improved  wettability  through  improved 
surface  energy  after  CNT  growth  [13].  Longer  growth  times  in  the  CVD  quartz  reactor 
improved  coverage  and  also  improved  wetting  to  the  point  where  the  PMMA 
became  fully  dispersed  on  the  surface.  Investigation  also  highlighted  the  kinetics  of 
infiltration;  the  resin  flowed  through  narrow  channels  between  the  CNT's,  with 
increasing  resin  viscosity  indicating  a  slower  diffusion  rate.  Here,  strong  bonding 
between  the  CNT  and  the  fibre  was  also  shown  through  a  fibre  fragmentation  test, 
reducing  fibre  fragmentation  lengths  indicative  of  improved  IFSS  [13]. 

CNT  growth  on  the  surface  of  preformed  fibres  in  a  composite  has  also  shown 
a  through  thickness  property  improvement;  something  which  is  desirable  for 
damage  tolerance  and  reductions  in  interlaminar  failure.  CNT  growth  on  a  3D  silicon 
carbide  weave  showed  interlaminar  fracture  toughness  improvements  of  348%  in 
mode  I  and  54%  in  mode  II  [34],  as  shown  via  illustration  in  Figure  2.  In  a  different 
experiment  on  a  unidirectional  carbon  fibre  composite,  flexural  strength  as  well  as 
modulus  were  improved  by  20%  and  28%  respectively  for  a  CNT  coated  preform 
[35],  Dependent  on  the  length  of  the  nano  growth,  they  could  be  considered  as  both 
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intra-laminar  and  inter-laminar  reinforcement  providing  resistance  to  delamination 
[34,  36],  Through  a  three  point  bend  test,  some  have  found  that  interlaminar  shear 
strength  is  increased  by  69%  after  growth  [37], 


Figure  4.2:  Illustration  of  'fuzzy'  fibres  through  radially  aligned  CNT  growth  on  a  woven  fibre  cloth 

[34] 

CNT's  have  also  been  grafted  to  the  surface  of  discontinuous  glass  fibres  and 
then  embedded  within  a  polypropylene  (PP)  resin.  It  was  found  that  the  in-plane 
tensile  strength  of  the  discontinuous  composite  was  improved  by  40%  compared 
with  a  bare  glass  fibre  PP  composite,  whilst  the  impact  strength  of  the  composite 
was  also  found  to  be  improved  by  24%  [38],  Post  failure  analysis  indicates  the  fibres 
pull  out  with  increased  adhesion  as  resin  still  remains  on  the  surface  (Figure  4.3). 
This  is  of  particular  interest  due  to  their  potential  application  in  ALM  through  the 
addition  of  short  fibres  to  the  feedstock. 


Figure  4.3:  SEM  images  of  a)  Smooth  surface  of  neat  glass,  b)  PP  sticking  to  the  fibre  surface  [38] 
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4.1.4  Zinc  oxide  nanorod  growth 

ZnO  nanorods  have  been  chosen  in  this  work  due  to  their  intrinsic  ductility 
and  higher  modulus  (29  GPa)  compared  with  typically  used  matrix  materials  (3-4 
GPa).  Compared  with  CNT's  (1,000  GPa),  the  lower  modulus  of  ZnO  can  still  provide 
load  transfer  between  discontinuous  fibres  while  allowing  ductile  failure.  ZnO 
nanorods  have  typically  been  grown  on  silica  substrates  previously  [18],  however, 
glass  has  a  similar  surface  chemistry  with  hydroxyl  groups  required  for  seeding  the 
substrate.  ZnO  nanorods  have  only  recently  been  grown  onto  curved  surfaces, 
primarily  optical  fibres  or  fibres  for  energy  harvesting,  shown  in  Figure  4.4  [22,  23], 
Current  uses  for  ZnO  nanorods  take  advantage  of  the  large  surface  area  of  nanorod 
arrays  as  compared  with  a  flat  substrate,  and  thus  efforts  have  been  focused  on 
tuning  growth  methods  for  the  densest  growth  morphologies.  If  the  dispersion  and 
alignment  of  the  rods  can  also  be  controlled,  IFSS  could  be  altered  and  controlled 
dependent  on  morphology. 


Figure  4.4:  SEM  image  of  ZnO  nanorod  coated  kevlar  fibre  for  energy  scavenging  applications  [22 ]. 

The  typical  two-stage  process  of  seeding  followed  by  growth  allows  for 
several  variable  factors  with  each  having  an  effect  on  morphology.  It  has  been 
shown  that  the  seeded  method  has  a  large  effect  on  morphology  [39,  40],  with 
sputtering  or  CVD  seeding  being  selected  for  more  control  over  nanorod  alignment 
as  well  as  producing  nanorods  with  higher  aspect  ratios  [22,  27,  41-44],  However, 
due  to  the  apparent  etching  effect  caused  by  the  CVD  or  sputtering  processes,  a  wet 
chemical  seeding  method  using  zinc  acetate  salts  has  been  used  here.  The  limitation 
of  quartz  reactor  size  or  sputtering  chamber  is  also  overcome  through  the  use  of  an 
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acetate  seed.  ZnO  is  expected  to  crystallise  through  the  hydrolysis  of  Zn  salts  in  a 
slightly  basic  solution,  and  hence  zinc  acetate  is  used  with  potassium  or  sodium 
hydroxide  in  order  for  crystallisation  to  occur  at  temperature  [18].  Ethanol  is 
typically  used  as  the  solvent  for  the  acetate  basic  solution  which  is  used  to  coat 
substrates  prior  to  growth.  Post  coating  of  the  substrate,  a  thermal  decomposition  is 
then  required  in  order  to  covalently  bond  the  ZnO  to  the  surface,  providing  a  glass 
surface  chemistry  change  shown  in  Figure  5.  Studies  have  shown  that  a  minimum 
temperature  for  thermal  decomposition  is  150°C,  with  higher  temperatures  around 
350°C  improving  perpendicular  alignment  through  a  reorientation  of  planes  [45], 
Seed  particles  in  ethanol  solution  vary  in  size  dependent  on  the  molarity  of  solution, 
temperature,  and  incubation  time  for  particles  to  form.  Typical  particle  sizes 
expected  for  a  5  mM  zinc  acetate  solution  are  5-15  nm  and  TEM  results  have  shown 
that  the  c-axis  of  the  crystal  wurtzite  structure  after  decomposition  on  to  a  substrate 
surface  (Figure  6)  are  oriented  perpendicular  to  the  surface  for  aligned  growth  [46], 
In  order  to  make  sure  a  substrate  is  fully  coated  with  zinc  acetate  seeds,  substrates 
are  typically  dipped  in  the  solution  several  times  followed  by  drying,  or  drop  cast  and 
spin  coated  onto  flat  substrates  several  times  [20,  23,  28,  39,  47];  with  the  latter  not 
a  feasible  option  for  cylindrical  fibres.  This  is  then  followed  by  a  high  temperature 
bake  to  encourage  particle-substrate  bonding. 

Once  seeded,  a  growth  method  can  be  selected  which  will  again  have  an 
influence  on  morphology.  More  control  is  seen  when  using  CVD  as  a  growth  method, 
allowing  for  much  higher  aspect  ratios,  with  rods  typically  much  longer  (Figure  7a), 
up  to  120pm  [48],  Flydrothermal  growth  on  the  other  hand  provides  the  prospect  of 
low  temperature  wet  chemical  synthesis  at  low  cost.  This  is  the  most  common 
method  in  literature  (Figure  7b)  [18],  and  unlike  alkali  mediated  growth  for  the  seed 
solution,  H20  is  the  source  of  oxygen  ions  required  for  ZnO  synthesis.  Low 
temperature  growth  and  seeding  allows  for  several  substrates  to  be  utilised  for 
growth  such  as  silicon  wafers,  glass,  thermoplastics  and  organic  fibres  [18].  An 
equimolar  solution  of  zinc  nitrate  hexamethylenetetramine  (HMTA)  is  used  in  order 
for  growth,  with  varying  molarity  across  literature  from  0.5  M  [47]  to  4  mM  [23],  In 
the  reaction,  zinc  nitrate  supplies  Zn2+  ions  whilst  HMTA,  a  non-ionic  cyclic  tertiary 
amine,  readily  hydrolyses  in  water  to  form  OH-  ions  to  form  a  basic  solution,  finally 
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combining  to  produce  ZnO  and  H20  after  several  intermediate  reactions  (Figure  8). 
The  exact  function  of  HMTA  is  still  relatively  unknown  [18],  however,  the  hydrolysis 
is  thought  to  occur  gradually,  causing  a  slow  release  of  OH-  ions  from  water,  with  a 
large  release  causing  ZnO  to  precipitate  out  too  quickly  causing  non-oriented 
growth.  Before  its  decomposition  to  produce  ions,  HMTA  attaches  to  the  non-polar 
surfaces  of  already  formed  nanorods,  promoting  longitudinal  growth.  Several  further 
structures  have  been  observed  through  seeding  after  initial  growth,  followed  by 
growth  for  a  second  time;  providing  secondary  smaller  rod  growth  perpendicular  to 
the  initial  rod  growth  [49],  Patterned  growth  can  also  be  observed  through  targeted 
seeding  of  the  surface  or  suppression  of  seed  nucleation  sites  [41].  This  report 
focusses  on  ZnO  nanorod  growth  on  glass  fibres  via  acetate  seeding  and 
hydrothermal  growth;  with  the  aim  of  improving  in-plane  properties  of 
discontinuous  composites  whilst  allowing  for  ductile  failure. 
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Figure  4.5:  Schematic  of  zinc  acetate  seeding  altering  the  surface  chemistry  of  glass. 
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Figure  4.6:  Wurtzite  crystal  structure  of  ZnO  nanorods. 
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Figure  4.7:  SEM  images  ofZnO  nanorods  a)  CVD  seed  and  growth  [48]  b)  acetate  seed  and 

hydrothermal  growth  [28]. 


HMTA  +  6H;0  <— >  4 NTH,  +  6HCHO 
NH3  +  H;0  <— >  NHj*  +  OH 
Zn:*+4NHj  »  [Zr^NH,)*)2* 


Zrr*  +  20H~  <— >  Zn(OH)2 


Zn(OH);  <— ►  Zn0  +  H;0 


Figure  4.8:  Intermediate  reactions  for  the  production  ofZnO  nanorods  during  hydrothermal  growth 

[18]. 

4.2  Experimental 

4.2.1  Materials  used 

Glass  microscope  slides  0.8-1  mm  thick  were  used  as  the  substrate  for 
morphology  control  on  a  flat  surface.  Discontinuous  short  chopped  glass  fibres  3  mm 
in  length  and  approximately  12  pm  in  diameter  from  vetrotextiles,  sized  with  a 
starch  based  oil  sizing  (T30),  were  used  for  batch  growth.  Continuous  glass  fibres 
were  produced  from  Duran  borosilicate  glass  with  a  diameter  of  approximately  32 
pm  with  no  sizing  agent  for  ZnO  growth.  Seed  solution  was  produced  with  zinc 
acetate  dihydrate  (Zn(CH3C00)2.2H20,  99%  purity.  Sigma  Aldrich),  ethanol 
(CH3CH20H,  absolute.  Sigma  Aldrich),  and  (KOH,  99%  purity.  Sigma  Aldrich).  For 
growth,  a  mix  of  zinc  nitrate  hexahydrate  (Zn(N03)2.6H20,  98%  purity.  Sigma 
Aldrich)  and  HMTA  (C6H12N4,  >99.0%  purity.  Sigma  Aldrich)  was  used  in  deionised 
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water  (MilliQ,  18.2MQ  resistivity,  <  3ppb  TOC).  Post  growth  a  tetraethyl  orthosilicate 
(TEOS)  (Si(OC2H5)4,  98%  purity.  Sigma  Aldrich)  coating  was  applied. 

4.2.2  Growth  method  ofZnO  on  flat  glass  surface 

Prior  to  growth  all  glass  slides  were  cleaned  in  69%  w/v  nitric  acid  (Sigma 
Aldrich).  Seeding  and  growth  was  performed  on  the  glass  microscope  slides  with  an 
approximate  area  of  1  cm2  post  nitric  acid  cleaning.  35.7  mg  of  Zn(CH3C00)2.2H20 
was  heated  for  lh  at  50°C  under  vigorous  stirring  in  40  mL  of  ethanol  (4  mM 
solution),  whilst  separately  12.9  mg  KOH  was  added  to  40  mL  of  ethanol  (4  mM 
solution)  and  heated  for  lh  at  50°C  under  vigorous  stirring.  Both  solutions  were  left 
to  cool  in  air  for  approximately  10  min  followed  by  mixing  through  drop  wise 
addition  under  vigorous  stirring.  The  4  mM  combined  solution  was  heated  at  70°C 
for  2  h  under  vigorous  stirring  and  then  bottled.  Seeding  was  performed  through 
dipping  the  cleaned  glass  slides  into  the  solution  for  15  min  at  a  time,  and  then 
placing  them  on  a  hot  plate  at  150°C  for  15  min  after  removal  from  the  solution.  This 
process  was  repeated  1-5  times  in  order  to  investigate  the  dependency  of 
morphology  on  the  density  of  the  seed.  After  coating  the  substrates,  they  were 
baked  at  310°C  on  a  hot  plate  for  2  h.  A  growth  solution  was  then  prepared  in  125 
mL  of  deionised  MilliQ  water  with  185.9  mg  of  Zn(N03)2.6H20  and  87.6  mg  of  HMTA 
(5  mM  solution).  Before  growth,  the  solution  was  stirred  vigorously  for  30  min. 
Substrates  were  then  placed  in  the  growth  solution  at  93°C  under  reflux  for  8-36  h 
(SI  Figure  1),  with  growth  time  varied  to  investigate  its  effect  on  morphology.  After 
growth,  substrates  were  removed  from  the  solution  and  rinsed  with  deionised  MilliQ 
water  to  remove  any  salts  and  heated  for  15min  at  150°C  to  remove  any  remaining 
water.  Finally  a  TEOS  coating  was  applied  through  dipping  the  substrate  in  a  TEOS 
solution  for  2-3  min  followed  by  drying  on  a  hot  plate  at  150°C. 

4.2.3  Growth  method  ofZnO  on  short  e-glass  fibre 

Short  e-glass  fibres  3  mm  in  length  initially  had  a  starch  based  oil  sizing  which 
was  removed  through  sonification  in  ethanol  for  1  h,  followed  by  three  15  min  UV 
ozone  treatments.  This  was  done  twice  to  ensure  polymer  removal.  This  opened  up 
hydroxyl  groups  on  the  glass  surface  required  for  bonding,  whilst  also  separating  the 
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fibres  after  they  were  bundled  together  post  sizing  (SI  Figure  2).  An  additional  step 
required  to  remove  the  fibres  from  solution  was  filtration  under  vacuum,  which 
occurred  after  sonification  (SI  Figure  3),  after  each  seeding  step,  and  after  growth. 
Creation  of  the  seed  and  growth  solution  remained  the  same  as  in  Section  2.2. 
Seeding  of  the  fibres  remained  constant  at  3  x  15  min  dipping  and  heating  cycles, 
followed  by  a  310°C  bake.  Growth  solution  concentration  was  varied  from  5  to  10 
mM  to  account  for  a  large  increase  in  growth  surface  area.  0.05  g  of  fibres  were 
batch  processed  at  a  time  in  order  to  control  surface  area  and  dispersion.  After 
growth,  fibres  were  removed  from  the  solution  and  rinsed  with  deionised  MilliQ 
water  to  remove  any  salts  and  then  heated  for  15  min  at  150°C  to  remove  any 
remaining  water.  Finally  a  TEOS  coating  was  applied  through  dipping  the  fibres  in  a 
TEOS  solution  for  2-3  min  followed  by  drying  on  a  hot  plate  at  150°C. 

4.2.4  Growth  method  ofZnO  on  continuous  fibre 

An  individual  borosilicate  glass  fibre  6cm  in  length  was  seeded  and  grown  on 
through  being  placed  in  a  mounting  fixture  and  clamped  at  either  end  (SI  Figure  4.4). 
Creation  of  the  seed  and  growth  solution  remained  the  same  as  in  Section  2.2.  Once 
in  the  mounting  fixture  it  remained  there  until  each  step  of  seeding  and  growth  had 
been  completed.  This  meant  that  growth  could  not  be  performed  in  a  round  bottom 
flask  under  reflux,  and  hence  had  to  be  performed  in  a  beaker  (SI  Figure  4.5).  A 
growth  solution  concentration  of  5  mM  was  used  for  8  h  of  growth.  Deionised  water 
was  added  to  the  solution  on  the  hour  in  order  to  maintain  growth  solution 
concentration.  After  growth,  the  fibril  was  rinsed  with  deionised  MilliQ  water  to 
remove  any  salts  and  heated  for  15  min  at  150°C  to  remove  any  remaining  water. 
Finally  a  TEOS  is  applied  by  dipping  the  fibre  in  a  TEOS  solution  for  2-3  min  followed 
by  drying  on  a  hot  plate  at  150°C. 

4.2.5  Individual  fibre  tensile  tests 

Single  fibre  tensile  tests  were  performed  at  room  temperature  in  accordance 
to  ASTM  C1557  [50],  with  a  gauge  length  of  10  mm  and  30  mm  (SI  Figures  4.6  and 
4.7).  Borosilicate  continuous  glass  fibres  were  tested  before  and  after  ZnO  nanorod 
growth  at  a  crosshead  speed  of  35  mm/min.  A  10  N  load  cell  was  used  on  a  1  kN 
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Instron  tensile  testing  machine,  whilst  a  Letoxit  epoxy  adhesive  (SI  Figure  4.6)  was 
used  to  attach  the  fibril  to  a  supporting  card.  Four  tests  were  performed  at  a  gauge 
length  of  30  mm  for  bare  glass  fibre,  eight  tests  at  10  mm  of  bare  glass  fibre  (SI 
Figure  4.8),  and  finally  four  tests  of  ZnO  nanorod  coated  fibres  following  treatment 
with  TEOS. 

4.2.6  Atomic  force  microscopy  measurements 

AFM  measurements  of  lateral  and  normal  force  were  taken  on  a  lx  seeded 
and  3x  seeded  ZnO  nanorod  coated  surface,  before  and  after  a  coating  of  TEOS.  It  is 
thought  that  this  would  be  indicative  of  the  strength  of  attachment  of  the  nanorods 
to  the  substrate  surface,  whilst  also  providing  some  insight  into  interfacial  shear  of 
different  nanorod  morphologies.  Tipless  Mikromasch  non-coated  CSC37  series  AFM 
probes  were  used  with  an  attached  glass  fibre,  both  with  and  without  ZnO  nanorod 
growth.  Lever  B  (SI  Figure  4.9)  with  a  force  constant  of  0. 1-0.6  N/m  was  selected. 
Tips  were  functionalised  through  particle  attachment  via  micromanipulators  using 
epoxy  glue  followed  by  curing.  Particles  attached  to  tipless  cantilevers  are  typically 
spherical  [51],  and  hence  this  method  will  provide  novel  insight  into  shape  effects  of 
the  cantilever  tip.  Normal  and  torsional  spring  constants  of  the  cantilevers  were 
determined  by  thermal  vibration  methods  [52,  53],  However,  this  method  showed 
significantly  lower  flexural  stiffness  than  expected  (SI  Table  2),  therefore  0.3N/m  was 
assumed  as  stated  as  the  average  by  the  Mikromasch  manufacturer  (SI  Figure  4.9). 
Tips  were  imaged  using  a  JEOL  JSM  6330F  for  high  resolution  images,  shown  in 
Figure  4.10.  Normal  force  measurements  were  performed  at  3-5  locations  across 
each  surface.  Lateral  force  measurements  (friction)  were  recorded  over  a  scan  range 
of  5  pm  at  1Hz,  providing  a  shear  velocity  of  10  pms-1,  with  variable  normal  force. 
Height  required  for  calculating  torque  was  assumed  as  particle  diameter  plus  half  the 
thickness  of  the  cantilever  (1  pm).  Higher  velocities  have  been  shown  to  cause  a 
reduction  in  friction  coefficient  [54],  likely  due  to  skipping  of  the  tip  over  the  rough 
surface.  Prior  to  use,  substrates  as  well  as  tips  were  cleaned  via  UV  ozone  treatment 
for  15  min. 
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Figure  4.9:  SEM  images  of  morphology  control  ofZnO  nanorods  through  altering  the  amount  of 
seeding:  a-c  show  1  repetition  of  the  seeding  process,  d-f  show  3  repetitions  of  the  seeding  process, 
g-i  show  5  repetitions  of  the  seeding  process  outlined  in  section  2.2.  Magnifications  remain  constant 

across  the  horizontal  axis. 
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Figure  4.10:  SEM  images  of:  a)  Glass  fibre  functionalised  tip  b)  Glass  fibre  with  ZnO  nanorod  growth 

functionalised  tip. 


4.3  Results  and  discussion 

4.3.1  Controlling  the  morphology  of  ZnO  growth  on  aflat  surface 

The  number  of  repetitive  seeding  cycles  on  a  1  cm2  substrate  was  varied  1-5 
times  followed  by  a  310°C  bake,  whilst  growth  conditions  initially  remained  constant 
with  a  5  mM  HMTA  and  zinc  nitrate  solution  at  93°C  for  8  h.  As  seen  in  Figure  4.9, 
increasing  the  seed  density  increased  nanorod  density  whilst  also  improving 
alignment.  It  is  also  seen  that  at  lower  seed  densities  nanorods  were  slightly  larger  at 
approximately  100  nm  in  diameter  compared  with  50  nm  for  the  3  and  5  times 
seeded  substrates.  This  is  thought  to  be  due  to  the  reduced  nucleation  sites  for 
growth  causing  increased  deposition  on  the  remaining  growth  sites.  Once  the  seed 
layer  reaches  a  certain  thickness,  rod  growth  impedes  adjacent  rods  from  forming  at 
an  angle  to  the  surface.  The  highest  density  is  likely  to  be  undesirable  for  wetting,  so 
a  3x  seeding  approach  was  selected  for  fibre  growth.  The  length  of  nanorods  was 
also  imaged  through  cracking  the  substrate  in  half,  as  shown  in  Figure  4.11,  to  be 
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approximately  1  pm  in  length.  It  was  found  that  varying  growth  time  at  8  h,  22  h  and 
38  h  had  no  effect  on  morphology  and  did  not  increase  nanorod  length;  indicating 
that  the  nutrients  in  solution  were  depleted  after  8  h  of  growth.  It  is  thought  that 
length  can  be  extended  through  increasing  concentration,  or  the  addition  of  further 
reactants  for  growth  after  a  certain  time  period. 


Figure  4.11:  SEM  image  of3x  seeded  substrate  grown  for  38  h  in  5  mM  growth  concentration. 

It  was  also  noted  that  morphology  varied  from  site  to  site  on  a  particular 
substrate,  as  shown  in  Figure  4.12.  It  can  be  seen  that  regions  of  higher  alignment 
are  well  pronounced  as  are  regions  with  lesser  alignment;  however,  with  increased 
seeding  the  regions  of  poorer  alignment  tend  to  diminish.  Variation  of  substrate 
surface  area  was  also  investigated,  indicating  a  morphology  dependence  that  is  so 
far  not  reported  in  literature.  Larger  rods  were  found  to  be  present  at  10  pm  in 
length  after  halving  the  substrate  surface  area;  however  they  were  not 
perpendicular  to  the  surface  and  not  clearly  bonded  to  the  surface  (SI  Figure  4.10). 

X-Ray  Diffraction  (XRD)  is  usually  performed  at  this  point  to  confirm  the 
degree  of  orientation  as  well  as  foreign  particle  content  [28];  however,  due  to  time 
constraints  this  was  not  performed.  A  high  resolution  image  clearly  shows  the 
hexagonal  structure  of  the  ZnO  nanorod  (Figure  4.13);  typical  of  the  wurtzite  crystal 
structure. 
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Figure  12:  SEM  image  of  4  times  seeded  sample,  with  lines  indicating  regions  of  different  alignment. 


Figure  13:  SEM  image  ofZnO  hexagonal  crystal  wurtzite  structure. 

4.3.2  ZnO  nanorods  on  glass  fibres 

After  initial  flat  substrate  testing,  a  cycle  of  3  seeding  steps  was  selected  for 
growth  on  short  glass  fibres  3  mm  in  length.  As  with  flat  glass  substrate  cleaning, 
concentrated  nitric  acid  cleaning  was  initially  employed  for  removal  of  the  starch 
based  polymer  sizing;  however,  it  was  found  to  damage  the  surface  of  the  glass  and 
caused  loss  of  nanorods  after  growth  (Figure  4.14).  During  this  stage  it  was  also 
noted  that  the  initial  5  mM  concentration  of  HMTA  and  zinc  nitrate  hexahydrate  for 
growth  was  insufficient  for  the  surface  area  requiring  growth  (SI  Table  1);  0.05  g  of 
fibres  has  an  approximately  20  times  increase  in  surface  area  compared  with 
previous  flat  substrate  growth.  The  growth  solution  was  therefore  increased  in 
concentration  to  10  mM,  whilst  sonification  in  ethanol  followed  by  UV  treatment 
was  selected  for  removal  of  the  polymer. 
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Good  coverage  was  noted  on  glass  fibres  after  an  increase  in  growth 
concentration  as  well  as  swapping  cleaning  methods  (Figure  4.15).  However,  fibres 
began  to  agglomerate  in  the  growth  solution  and  bundles  of  glass  fibres  formed; 
even  whilst  under  vigorous  stirring  (SI  Figure  4.11).  Central  areas  of  the  bundles 
show  poor  nanorod  growth,  while  outer  areas  exposed  to  nutrients  received 
excellent  coverage  and  growth.  Approximate  ZnO  nanorod  aspect  ratio  remained 
similar  to  that  on  flat  glass  substrates  and  perpendicular  to  the  surface  (SI  Figure 
4.12).  It  is  also  seen  from  Figure  4.15  that  there  is  an  increase  in  the  deposition  of 
the  nanorods  on  the  surface  that  have  formed  in  solution.  Flat  glass  slides  were 
easily  cleaned  with  MilliQ  water;  however,  filtration  under  vacuum  to  remove  fibres 
from  solution  also  removes  longer  nanorods  that  have  grown  in  solution.  This  makes 
cleaning  of  the  fibre  bundles  more  difficult;  however,  rinsing  was  still  conducted 
several  times  whilst  fibres  remained  on  the  fine  glass  Buchner  funnel.  It  is  very 
difficult  to  avoid  bundling  during  growth  as  open  hydroxyl  surface  groups  are 
required  for  ZnO  nanorod  surface  bonding,  and  also  cause  hydrogen  bonding  of  the 
glass  fibres. 


Figure  4.14:  SEM  image  of  the  loss  of  ZnO  nanorods  from  the  glass  fibre  surface  due  to  concentrated 

acid  cleaning. 
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Figure  4.15:  SEM  image  of  3x  seeded  short  glass  fibre,  10  mM  growth  concentration  in  a  batch 

process. 

In  order  to  counter  the  agglomeration  problem  as  well  as  any  unforeseen 
consequences  of  polymer  sizing  removal,  continuous  borosilicate  glass  fibres  6  cm  in 
length  were  produced  without  any  sizing.  A  larger  diameter  was  produced  of 
approximately  33.28  pm  in  diameter  (SI  Figure  4.13),  where  before  and  after  growth 
can  be  seen  in  Figure  4.16.  Distribution  is  noted  to  be  sparser  than  previously  grown 
rods  on  both  a  flat  glass  surface  and  short  glass  fibres  (SI  Figure  4.17).  This  is  thought 
to  be  due  to  the  poor  temperature  seeding;  the  rig  setup  does  not  allow  an  oven 
bake  and  cannot  be  placed  in  direct  contact  with  the  hot  plate  surface.  Therefore, 
although  the  hot  plate  temperature  was  set  at  a  maximum,  the  fibre  temperature 
was  considerably  lower.  The  quality  of  the  seed  has  been  previously  shown  to  be 
dependent  on  temperature  [45],  and  hence  a  custom-designed  holding  device  will  be 
required  to  further  the  potential  of  this  method  that  allows  temperatures  up  to 
350°C  within  a  furnace.  Due  to  the  fewer  nucleation  sites,  nanorods  have  grown  with 
an  approximate  aspect  ratio  of  10;  considerably  less  than  the  aspect  ratio  of  20  seen 
on  flat  glass  substrates  and  short  glass  fibres.  Contaminants  can  also  be  seen  on  the 
fibre  surface,  which  is  thought  to  be  due  to  the  oxidation  of  the  clamping 
equipment. 
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Figure  4.16:  a)  SEM  image  of  glass  fibre  prior  to  ZnO  nanorod  growth  indicating  a  fibre  diameter  of 

32.81pm,  b)  after  ZnO  growth. 


Figure  4.17:  SEM  image  of  continuous  fibre  ZnO  nanorod  growth,  indicating  sparser  arrangement  of 

nanorods. 

In  order  to  ensure  the  ZnO  nanorods  are  securely  seated  on  the  fibre  surface, 
a  TEOS  coating  is  applied.  TEOS,  often  used  as  a  cross  linking  agent  or  to  make  pure 
silica  [55]  (Figure  4.18),  cross  links  on  the  surface  when  combined  with  moisture  in 
the  atmosphere.  This  in  effect  embeds  the  roots  of  the  ZnO  fibre  to  the  glass  surface 
in  silica  [22]  (Figure  4.19).  Through  increasing  surface  roughness  and  contact  area, 
nanorods  or  nanotubes  are  expected  to  draw  the  non-viscous  TEOS  to  the  roots 
through  capillary  action  [13,  22,  29],  A  brief  investigation  of  SEM  before  and  after 
TEOS  coating  revealed  no  difference  in  imaging,  ensuring  that  TEOS  had  not  been  left 
on  the  surface  covering  nanorods  (SI  Figure  4.14-4.15). 
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Figure  4.18:  TEOS  cross  linking  in  the  presence  of  H20. 
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Figure  4.19:  Illustration  of  TEOS  coated  fibre  cross  section;  not  to  scale. 


4.3.3  Individual  fibre  tensile  testing 

Individual  fibre  tensile  tests  were  performed  on  the  produced  continuous 
borosilicate  glass  fibre,  approximately  33.281  pm  in  diameter  (SI  Figure  13),  before 
and  after  ZnO  nanorod  growth.  This  was  performed  to  understand  if  the  attachment 
method  had  any  unexpected  detrimental  effects  on  tensile  strength  similar  to  that  of 
the  CVD  approach  seen  for  CNT  growth.  An  initial  gauge  length  of  30  mm  was 
chosen,  followed  by  10  mm  in  order  to  reduce  the  probability  of  premature  failure. 
The  results  of  tensile  strength  are  given  in  Figure  20.  The  lowest  tensile  strength 
average  of  278  MPa  was  seen  at  a  gauge  length  of  30  mm,  whilst  reducing  gauge 
length  to  10  mm  increased  the  average  tensile  strength  to  336  MPa.  As  seen  from 
the  SEM  in  Figure  16(a),  there  are  several  defects  in  the  glass  produced,  having  a 
pronounced  effect  on  the  fibre  tensile  strength  [14].  By  reducing  gauge  length  by  a 
third  the  probability  of  inclusion  of  a  limiting  strength  defect  in  the  glass  is  also 
reduced  by  a  third.  At  a  gauge  length  of  10  mm,  the  average  strength  found  differed 
by  only  4%  from  literature  where  identical  glass  and  production  procedures  were 
used  [56],  Tensile  strength  following  ZnO  growth  and  a  TEOS  coating  showed  an 
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average  increase  to  373  MPa;  clearly  indicating  that  the  growth  procedure  is  not 
detrimental  to  the  tensile  properties  of  the  borosilicate  glass  but  in  fact  provides  an 
approximate  11%  increase  in  strength.  This  is  thought  to  be  due  to  the  TEOS  coating 
following  growth;  where  the  TEOS  coating  'heals'  some  of  the  imperfections  shown 
in  Figure  16(a),  reducing  stress  concentration  factors  imparted  on  the  fibre  from  the 
defects.  Also  expected  is  a  slight  cross  sectional  area  increase  (Figure  19);  however  a 
3.3%  increase  in  radius  relating  to  an  11%  increase  in  area  seems  unlikely.  In  order  to 
confirm  this,  a  further  comparative  study  would  be  required  on  the  tensile  strength 
of  bare  glass  fibres  with  just  a  TEOS  coating.  Young's  modulus  was  also  calculated  to 
be  65  GPa  after  ZnO  growth,  slightly  higher  than  anticipated  with  Duran  Borosilicate 
glass  expected  to  have  a  modulus  of  63GPa.  Raw  data  is  given  in  SI  Figures  14-16. 


Tensile  Strength  Comparison 


Figure  4.20:  Comparative  single  fibre  tensile  strengths  at  30  mm  gauge  length,  10  mm  gauge  length 
and  10  mm  gauge  length  with  ZnO  nanorod  growth  and  a  TEOS  coating. 

4.3.4  AFM  measurements 

Normal  force  data  was  taken  to  ensure  the  cantilever  was  operating  within  its 
linear  range  of  the  photodiode  and  inspect  any  adhesion  forces  present.  This  data 
also  provided  the  cantilever  sensitivity  (nm/V)  and  the  zero  load  offset  for 
deflection;  calculated  each  time  the  AFM  tip  probe  was  re-engaged  to  the  substrate 
surface.  Glass  against  a  bare  glass  fibre  was  analysed  first  for  a  base  comparison, 
indicating  good  separation  of  trace  and  retrace  curves  (Figure  4.21).  Roughness  of 
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the  surface  is  shown  through  peaks  and  troughs;  with  the  trace  and  retrace 
appearing  similar  (i.e.  of  mirror  image  of  each  other),  suggesting  little  skipping  of  the 
probe  due  to  high  velocities.  A  coefficient  of  friction  was  calculated  according  to 
Amontons'  laws  of  friction  [57],  taking  an  average  lateral  force  from  the  trace  and 
retrace  from  16  loading  conditions  from  0  to  95  nN.  Friction  can  then  be  plotted 
against  normal  load  and  a  linear  best  fit  applied  to  provide  the  friction  coefficient 
(Figure  4.22).  Figure  4.22  shows  a  difference  between  loading  and  unloading 
averages  which  tend  to  increase  at  lower  normal  force  loading.  This  may  be  due  to  a 
smaller  normal  force  leading  to  a  slower  reaction  and/or  poor  surface  contact 
following  asperities  on  the  surface.  A  friction  coefficient  of  3.24  was  calculated, 
higher  than  previously  published  results  of  AFM  tip's  against  pure  silicon  of  1.07  [54], 
It  is  thought  that  this  is  due  to  errors  in  the  spring  constant  stiffness  and  torsional 
calculation  (which  will  be  checked)  however,  results  compared  here  will  be  relative 
and  thus  internally  consistent. 

The  3x  seeded  substrates  were  then  investigated  prior  to  and  after  a  TEOS 
coating  in  order  to  compare  adhesion.  Particle  contact  area  estimation  is  difficult 
due  to  the  angle  of  the  tip,  its  shape,  and  approximate  length;  and  hence  only  an 
estimate  of  contact  pressure  is  given  here.  Assuming  a  contact  area  of  a  6th  of  the 
projected  area  of  the  cylinder,  a  contact  pressure  of  approximately  3.4  kPa  is  applied 
at  a  maximum  load  of  95  nN.  This  is  very  small  compared  with  stress  expected  to  be 
transferred  through  IFSS  in  a  composite;  however  true  contact  area  is  expected  to  be 
lower  due  ZnO  nanorod  asperities.  This  causes  regions  of  higher  local  stress,  whilst 
the  differences  in  the  trace  and  retrace  have  potential  to  indicate  a  possible 
direction  dependence  of  friction  due  to  anisotropic  alignment  of  rods.  It  is  also 
possible  that  a  difference  in  trace  and  retrace  is  due  to  ZnO  nanorod  cleavage  from 
the  surface.  Raw  traces  at  95  nN  of  applied  load  of  3x  seed  without  and  with  TEOS 
can  be  seen  in  Figure  4.23  and  Figure  4.24  respectively.  Both  traces  indicate  a  much 
rougher  surface  with  a  large  increase  in  the  amplitude  of  peaks  and  toughs  as 
expected  after  ZnO  nanorod  growth  [54],  Large  amplitude  peaks  of  up  to  8V  were 
often  seen  in  both  the  trace  and  retrace  largely  due  to  rods  that  were  not  aligned 
perpendicular  to  the  surface.  An  increase  in  the  asymmetry  of  peaks  of  amplitudes  at 
approximately  3 V  was  often  seen  in  3x  seeded  growth  without  TEOS,  perhaps 
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indicating  damage.  Furthermore,  friction  plots  of  the  substrate  without  TEOS  tend  to 
show  a  drop  in  friction  on  unloading  providing  hysteresis,  again  indicative  of  damage 
(Figure  4.25).  Following  a  TEOS  coating,  hysteresis  is  generally  eliminated  (Figure 
4.26).  This  is  shown  at  several  locations  across  the  surface  (SI  Figure  4.19).  The 
coefficient  of  friction  on  average  was  found  to  be  6.84  following  a  TEOS  coating; 
directly  comparable  to  the  uncoated  ZnO  nanorod  growth  friction  coefficient  of  6.97. 
Following  these  results  all  substrates  were  coated  with  TEOS  for  the  remainder  of 
testing. 

The  inter-digitation  of  nanorod  coated  surfaces  was  also  of  interest.  A 
nanorod  coated  fibre  with  3x  seeded  growth  was  used  as  a  colloidal  particle  to 
perform  a  similar  test  on  the  3x  seeded  substrate,  both  coated  in  TEOS.  Shear  traces 
show  good  separation  of  surfaces  indicating  friction  was  substantial  (Figure  4.27). 
Peaks  and  troughs  are  of  a  similar  frequency  and  amplitude  as  in  previous  shear 
traces.  The  friction  plots  show  little  to  no  hysteresis  indicating  rods  are 
perpendicular  to  the  substrate  surface  and  no  damage  is  observed  on  the  retrace. 
Friction  coefficient  is  calculated  as  5.49,  slightly  lower  than  previous  measurements 
of  ZnO  nanorod  growth  against  a  bare  glass  fibre.  This  is  thought  to  be  due  to  the 
inter-digitation  of  nanorods  causing  a  reduction  in  nanorod  bending,  in  turn  reducing 
friction.  A  smoother  surface  on  the  other  hand  has  contact  with  more  rods  causing 
increased  rod  bending.  A  reduction  in  the  coefficient  of  friction  between  nanorod 
coated  surfaces  has  been  seen  previously  supporting  this  result  [54], 

The  lx  seeded  ZnO  growth  substrates  with  lesser  alignment  were  also  used  in 
AFM  with  the  hope  of  identifying  a  friction  dependence  on  nanorod  morphology. 
Flowever,  results  showed  much  lower  friction  coefficients,  thought  to  be  due  to 
contact  with  larger  rods  formed  in  solution  connected  to  the  surface  (SI  Figure  4.20). 
It  is  thought  that  this  causes  1  pm  length  nanorods  to  be  out  of  contact,  therefore 
not  contributing  to  the  frictional  resistance  seen.  Thorough  analysis  has  however  yet 
to  be  performed  on  this  data. 
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Figure  4.21:  Bare  glass  substrate  vs  bare  glass  probe  shear  trace  following  95  nN  of  force  assuming 


a  spring  constant  of  0.3  N/m 


Figure  4.22:  Frictional  plot  of  bare  glass  substrate  vs  bare  glass  probe  with  varying  normal  load  with 
a  linear  line  of  best  fit  (Amonton's  law),  suggesting  a  friction  coefficient  of  3.24. 


Figure  4.23:  3x  seeded  ZnO  vs  bare  glass  probe  shear  trace  following  95  nN  of  force  assuming  a 

spring  constant  of  0.3  N/m. 
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Figure  4.24:  3x  seeded  ZnO  vs  bare  glass  probe  shear  trace  following  95  nN  of  force  assuming  a 


spring  constant  of  0.3  N/m. 


Figure  4.25:  Frictional  plot  of  a  3x  seeded  ZnO  nanorod  substrate  vs  bare  glass  probe  with  varying 
normal  load  with  a  linear  line  of  best  fit  (Amontons'  law),  giving  a  friction  coefficient  of  6.97 


Figure  4.26:  Frictional  plot  of  a  3x  seeded  ZnO  nanorod  substrate  +  TEOS  vs  bare  glass  probe  with 
varying  normal  load  with  line  of  best  fit  (Amontons'  law),  with  a  friction  coefficient  of  6.84. 
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Figure  4.27:  3x  seeded  ZnO  +TEOS  vs  3x  seeded  ZnO  +TEOS  glass  fibre  probe  shear  trace  (95nN 

assuming  0.3  N/m ) 


Figure  28:  Frictional  plot  of  a  3x  seeded  ZnO  nanorod  substrate  +  TEOS  vs  3x  seeded  ZnO  nanorod 
substrate  +  TEOS  glass  fibre  probe  with  varying  normal  load  with  a  linear  line  of  best  fit  (Amontons' 

law),  giving  a  friction  coefficient  of  5.49. 

4.4  Concluding  Remarks 

Morphology  and  alignment  of  ZnO  nanorod  growth  through  hydrothermal 
methods  combined  with  an  acetate  seeding  approach  has  been  shown  to  be  highly 
dependent  on  the  seed  density.  Increasing  seed  density  improves  alignment 
perpendicular  to  the  surface  through  the  restriction  of  growth  on  non-aligned  axes. 
Batch  processed  ZnO  nanorod  growth  on  short  glass  fibres  has  been  shown  to  be 
successful;  however  only  in  areas  open  to  nutrient  supply  near  the  perimeter  of 
agglomerated  fibres.  Agglomeration  of  discontinuous  fibres  during  growth  limits  the 
potential  for  a  batch  process  to  be  utilised  for  composite  fibre  manufacture  with  ZnO 
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nanorod  attachments,  and  a  continuous  fibre  growth  method  where  fibres  can  be 
cut  after  growth  is  recommended.  The  method  of  cleaning  has  also  shown  to  have 
had  a  significant  effect  on  adhesion  to  the  surface;  with  concentrated  acid  damaging 
the  fibre  surface  causing  nanorod  loss.  Also  shown  but  not  fully  investigated  is  the 
effect  on  growth  morphology  of  variable  substrate  surface  area,  something  yet  to  be 
reported  in  current  literature.  Growth  on  a  continuous  glass  fibre  through  clamping 
at  either  end  is  shown  to  be  successful  with  extremely  good  coverage  of  nanorods. 
Poor  heating  of  the  individual  fibril  during  decomposition  of  the  zinc  acetate  seed  is 
thought  to  be  the  cause  of  sparser  rod  density,  whilst  surface  contaminants  and 
poorer  heating  during  growth  are  also  thought  to  contribute.  Tensile  strength  after 
ZnO  nanorod  growth  followed  by  a  TEOS  coating  was  found  to  increase,  with  the 
increase  likely  due  to  the  TEOS  coating  'healing'  surface  defects  and  effectively 
increasing  effective  cross  sectional  area.  AFM  measurements  have  shown  that  after 
a  TEOS  coating  there  is  reduced  lateral  force  hysteresis,  potentially  indicating  less 
surface  damage  from  the  AFM  probe.  Inter-digitation  of  nanorods  against  nanorods 
has  shown  a  reduction  in  friction  compared  with  nanorods  against  a  bare  glass 
surface,  thought  to  be  due  to  reduced  rod  bending  caused  by  an  overall  reduction  in 
contact  area.  The  overriding  observation  from  the  friction  measurements  is  that  the 
presence  of  nanorods  has  a  dramatic  effect  on  the  interfacial  shear,  which  can  be 
tuned  by  controlling  the  nanorod  morphologies.  This  paves  the  way  for  future 
exploitation  of  nanorod  coating  as  an  effective  strategy  for  controlling  interfacial  slip 
and  stress  in  composites  for  enhanced  ductility. 

There  are  several  areas  for  future  research  for  ZnO  nanorod  application  in 
composite  materials.  First  and  foremost  is  to  quantify  IFSS  of  a  fibre  post  surface 
nanorod  growth,  either  through  fibre  pull-out  or  a  fibre  fragmentation  test.  Further 
to  this,  potential  for  growth  on  an  individual  continuous  preform  is  possible 
providing  the  preform  has  no  sizing  agent.  Due  to  its  crystal  wurtzite  structure,  ZnO 
nanorods  have  piezoelectric  behaviour  and  hence  could  be  used  as  energy 
harvesters  and/or  as  strain  sensors.  Carbon  fibres  are  naturally  conducting  and  may 
provide  potential  to  carry  charge  produced  via  induced  strain;  providing  the  enticing 
prospect  of  smart  fibres.  Finally,  patterned  growth  through  different  seeding 
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techniques  is  possible,  which  may  provide  potential  for  IFSS  control,  controlling 
where  facture  may  initiate  and  propagate. 
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5.  Ultrasonic  alignment  of  micron  length  scale 
reinforcement 


5.1  Introduction 

To  this  day,  the  majority  of  ultrasonic  devices  have  been  developed  for 
applications  in  biological  and  life  sciences.  Consequently,  designs  are  generally 
tailored  to,  and  optimised  for  the  operation  in  an  aquaeous  environment.  In  this 
study,  the  primary  objective  is  to  develop  a  design  that  allows  for  the  rapid  and 
repeatable  fabrication  of  thin  layers  of  composite  material  at  the  micron  length 
scale;  a  secondary  motivation  is  founded  on  the  desire  to  maximise  the  range  of 
fibre  architectures  that  can  be  formed  using  a  particular  ultrasonic  device.  The  final 
design  is  to  be  devised  taking  a  counter-propagating  wave  approach,  as  this  method 
is  considered  to  be  the  least  sensitive  to  changes  in  the  fluid  chamber's  resonant 
frequency  when  large  quantities  of  particles  (at  the  micron  length  scale)  are  present. 
Moreover,  samples  are  to  be  sufficiently  large  to  allow  for  macroscopic  mechanical 
characterisation  to  be  carried  out. 

5.2  Device  design:  Two-element  array 

The  two-element  array  is  the  simplest  form  of  counter-propagating  wave 
device  and  was  used  extensively  in  this  study  for  the  manufacture  of  thin  composite 
samples  with  unidirectional  reinforcement.  The  device  is  shown  schematically  in 
Figure  5.1  below. 
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Figure  5.1:  Third  generation  two-element  array,  (a)  A  standing  wave  acoustic  field  is 
formed  inside  a  central  fluid  volume,  by  interference  of  two  counter-propagating 
traveling  waves.  The  PZT  transducers  are  submerged  in  water  for  cooling  purposes, 
and  held  in  place  by  small  compression  springs,  (b)  A  device  representation  by 
material  layers:  air  (1),  PMMA  (5  mm),  H20  (9.025  mm),  PZT  (0.975  mm),  PMMA 
(5mm),  resin  (30mm),  PMMA  (5mm),  PZT  (0.975mm),  H20  (9.025mm),  PMMA 
(5mm),  air  (1). 

As  shown  in  Figure  5.1,  two  0.975  mm  x  15  mm  x  2  mm  soft-doped  PZT 
(Noliac  group,  NCE51)  elements  were  placed  on  opposite  ends  of  the  device,  each 
separated  from  the  central  cavity  by  a  sacrificial  PMMA  boundary.  On  either  side  of 
the  device,  an  additional  water  filled  chamber  served  the  purpose  of  a  necessary 
heat  sink  at  high  driving  voltages.  The  PMMA  frame  was  mounted  on  a  glass 
substrate  using  adhesive  tape  (tesa  64621-00007-01).  For  both  acoustic  drivers  to  be 
more  easily  recycled,  a  small  compression  spring  secured  each  of  them  against  the 
adjoining  PMMA  boundary;  a  significantly  increased  composite  production  efficiency 
was  the  immediate  result.  With  the  central  cavity  now  measuring  30  mm  x  15  mm  x 
2  mm  (a  size  suitably  large  to  manufacture  test  specimens  for  mechanical 
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characterization),  the  driving  voltage  was  raised  to  80  Vpp.  An  operation  at  relatively 
low  voltages  is  desirable  primarily  to  minimise  the  effects  of  acoustic  streaming 
inside  the  central  cavity.  In  the  x-y  plane  the  sample  size  was  constrained  to  the 
dimensions  of  the  device  cavity,  the  thickness  (z-direction)  could  be  controlled  by 
varying  the  amount  of  resin  added  to  the  chamber. 

The  composite  assembly  time  is  critically  dependent  on  viscosity,  with  the 
fluid  drag  constituting  the  dominant  resistive  force  to  any  particle  motion.  A  low 
viscosity  resin  is  thus  generally  preferred.  For  short  glass  fibres  (r  <  25  pm)  dispersed 
in  water  (nwater  ~  0.8  x  10-3  Pa  s),  full  assembly  is  experimentally  achieved  in  less 
than  100  ms,  with  particles  traveling  to  pressure  nodes.  Given  a  linear  dependency 
of  the  drag  force  on  viscosity,  an  equivalent  assembly  time  of  950  ms  is  calculated 
for  fibres  suspended  in  the  resin. 

5.3  Materials 

A  blend  of  fibre  composite  material  used  was  based  on  the  aromatic 
dimethacrylate  monomer  2,2-bis[p-(2-hydroxy-3- 

methacryloxypropoxy)phenyl]propane  (BisGMA),  a  foundation  commonly  used  for 
dental  composite  applications.  To  lower  the  viscosity  of  the  bulky  BisGMA  monomer 
(from  r|BisGMA  ~  800  Pa  s  to  nreSin  ~  0.76  Pa  s  [24]),  it  was  mixed  with  the  comonomer 
triethylene  glycol  dimethacrylate  (TEGDMA)  in  a  ratio  of  1:1  by  weight.  0.2wt%  of 
camphorquinone  (CQ)  and  0.8wt%  of  Ethyl  4-  (dimethylamino)benzoate  (EDB)  were 
also  added  to  achieve  radical  photo-polymerisation  at  a  wavelength  of  460  nm.  Table 
5.1  illustrates  the  molecular  structures  of  each  of  these  components.  Complete 
ploymerisation  was  reached  following  exposure  to  a  900  mW  light-emitting  diode 
(LED)  for  approximately  40  s. 
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Table  5.1:  Molecular  structures  of  BisGMA,  TEGDMA,  EDB,  and  CQ. 


CH3  OH  OH  ch3 

2,2-bis(p-(2-hydroxy-3-methacryloxypropoxy)phenyl)propane  (BisGMA) 


CH3  O 

Triethylene  glycol  dimethacrylate  (TEGDMA) 


H,C 


Ethyl  4-(dimethylamino)benzoate  (EDB) 


Camphorquinone  (CQ) 


Commercially  available,  milled  glass  fibres  (Lanxess,  MF7904)  were  used  as 
reinforcement,  with  a  nominal  length  of  50  pm  and  approximately  14  pm  in 
diameter. 


5.4  Fabrication  of  epoxy  test  specimens 

The  fabrication  process  for  the  ultrasonically  assembled  samples,  is  shown 
diagrammatically  in  Figure  5.2.  Firstly,  some  dry  fibre  reinforcements  were 
distributed  on  the  substrate.  About  0.8  ml  of  resin  were  then  added,  before  the  two 
phases  were  mixed  to  give  an  even  dispersion  of  glass  particles  on  the  bottom 
surface.  An  ultrasonic  field  was  applied  to  drive  the  fibres  to  the  nodal  positions  of 
the  resulting  standing  wave.  Given  a  driving  frequency  of  v  =  2.01  MHz  and  assuming 
a  speed  of  sound,  cl  =  1480  m  s-1  in  the  fluid,  lines  of  fibres  are  predicted  to  form  at 
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a  separation  of  A/2  =  368  pm;  experimental  measures  find  this  separation  to  be  in 
the  region  of  350  pm  to  380  pm.  Prior  to  polymerisation  initiation,  reinforcement 
entities  were  allowed  approximately  60  s  to  arrange  and  orient  themselves  under 
the  influence  of  ultrasound.  Finally,  the  resin  was  exposed  to  a  blue  light  LED  and 
cured  in  40  s.  The  microstructure  of  the  resulting  samples  is  shown  in  Figures  5.3  (a) 
and  (e). 


Figure  5.2:  Sample  fabrication  process. 
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Figure  5.3:  Ultrasonic  assembly  of  glass  particles  into  lines  of  unidirectional 
reinforcement  in  (a)  epoxy  (c  =  1510  m  s'1  ±  2.8  %),  (b)  polyester  (c  =  1347  m  s'1  ±  2.5 
%),  (c)  silicone  (c=1070ms~1  ±4.4%),  (d)  agar  (c=1513ms~1  ±2.9%).  Examples  are 
shown  at  different  magnifications  and  for  varying  fibre  volume  fractions.  A  typical 
sample  cross-section  (e)  is  presented  in  the  form  of  a  x-ray  micro-tomography  image. 
Fibre  lines  extend  along  the  y-direction  and  are  separated  in  the  x-direction  by  a 
theoretical  distance  of  A/2  =  c/2f  where  c  is  the  speed  of  sound  in  the  matrix  material 
and  f  is  the  driving  frequency. 

Following  the  removal  of  the  acoustic  drivers  and  the  substrate,  two  identical 
13.6  mm  x  13.6  mm  test  pieces  were  cut  from  the  original  sample  using  a  high 
precision  carbon  dioxide  (C02)  laser.  Across  samples,  the  fibre  volume  fraction  was 
estimated  to  be  Vf  =  9.0  ±  2.0  %,  by  pCT  analysis;  for  each  pair  of  test  pieces,  Vf  was 
assumed  to  be  the  same.  By  means  of  the  above  procedure  a  true  comparison  of 
mechanical  test  results  could  be  established  for  each  pair  of  test  pieces.  An  analysis 
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between  samples  also  serves  to  quantify  general  manufacture  consistency,  to 
estimate  average  strength  and  stiffness  values,  and  to  determine  failure 
mechanisms.  The  slight  meniscus  formed  by  excess  resin  on  top  of  each  specimen 
was  removed  by  abrading  all  specimens  to  a  uniform  thickness  of  0.50  ±  0.05  mm.  To 
allow  for  better  gripping  of  the  test  pieces  during  mechanical  testing  (in  tension), 
end  tabs  -  as  depicted  in  Figure  5.4  -  were  attached  to  the  specimens.  A  set  each  of 
randomly  dispersed  fibre  samples  and  neat  resin  specimens  were  prepared  in  the 
same  manner  but  no  ultrasonic  energy  was  applied  in  these  cases. 


Figure  5.4:  Schematic  illustration  of  the  specimen  geometry  for  mechanical 
characterisation.  End  tabs  were  adhered  to  each  pair  of  test  specimens.  In  one  case 
end  tabs  were  attached  parallel  to  the  direction  of  fibre  reinforcement,  in  the  other 
transversely  oriented. 

5.5  Mechanical  Characterisation 

The  samples  were  mechanically  tested  in  tension  at  a  displacement 
controlled  crosshead  speed  of  0.10  mm  min-1.  The  results  are  summarised  in  the 
form  of  stress  vs.  strain  plots,  shown  in  Figure  5.5  with  the  average  failure  and 
stiffness  results  summarized  in  Table  5.2. 


R.  S.  Trask  @  Bristol,  ac.  uk 


67 


30/03/2014 


Distribution  A:  Approved  for  public  release;  distribution  is  unlimited. 


Additive  Layer  Manufacturing  Of  Biologically  Inspired  Short  Fibre  Reinforced  Composites 

Prepared  by  Dr  R.S.  Trask 


a-  b. 


Figure  5.5:  A  comparison  of  (a)  stiffness,  and  (b)  failure  stress  between  samples 
undergoing  tensile  loading  in  the  x-  and  in  the  y-direction.  For  both  cases,  the  mean 
and  standard  error  are  indicated  as  a  reference. 


Table  5.2:  Average  results  for  strain  to  failure,  stress  to  failure,  and  stiffness.  Values 
are  presented  for  four  different  sample  types,  and  the  associate  standard  error  is 

given  in  all  cases. 


Sample  type  Failure  strain,  e 

_ [%] 

No  fibres  4.55  ±  0.41 

Random  fibre  distribution  2.63  ±  0.17 

^-direction  2.04  ±  0.13 

y-direction  3.00  ±  0.15 


Failure  stress,  a 

[MPa] 

35.07  ±  1.24 
31.76  ±  1.85 
31.64  ±  1.61 
45.05  ±  1.97 


Stiffness,  E 

[MPa] 
11.51  ±  0.61 
13.33  ±  0.61 
16.36  ±  0.48 
17.66  ±  0.63 


An  average  stiffness  of  16.36  ±  0.48  MPa  was  recorded  for  transversely 
strained  samples;  in  comparison,  specimens  loaded  in  the  direction  of  fibre 
reinforcement  consistently  showed  superior  mechanical  performance  with  an 
average  stiffness  of  17.66  ±  0.63  MPa.  The  Young's  modulus  of  samples  without 
reinforcement  and  randomly  oriented  reinforcement  were  obtained  as  11.51  ± 
0.61MPa  and  13.33  ±  0.61MPa,  respectively.  It  is  thus  evident  that  improved 
structural  properties  and  anisotropy  have  been  achieved.  On  comparison  between 
samples  with  randomly  oriented  reinforcements  and  ultrasonically  assembled 
samples  strained  transversely,  it  should  be  noted  that  the  latter  appear  to  provide 
greater  stiffness.  It  is  believed  that  this  is  due  to  the  fact  that  within  the  assembled 
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lines,  fibres  are  in  close  proximity  allowing  interaction  and  load  transfer,  which 
results  in  increased  stiffness.  Within  the  random  samples  the  fibres  are  dispersed 
and  cannot  interact. 

As  is  clear  from  both  Figure  5.5  and  Table  5.2,  the  matrix  itself  exhibits  a 
more  ductile  mechanical  behaviour  with  strains  to  failure  of  up  to  5.55  ±  0.41  %;  on 
the  addition  of  fibre  reinforcement,  a  more  brittle  failure  mechanism  is  observed 
with  failure  generally  occurring  in  a  region  of  2  to  3  %  strain. 

For  an  estimation  of  the  expected  composite  stiffness,  Ey,  the  basic  rule  of 
mixtures  (assuming  continuous  fibre  reinforcement)  was  applied: 

Ey  =VfEf  +(l-Vf)Em. 

Flere,  Em,  and  Ef  are  the  Young's  modulus  of  the  matrix  and  the  fibres, 
respectively.  Across  all  samples,  an  average  Ey  =  18.97  ±  0.57  MPa  was  so  calculated, 
taking  Vf  =  9.0  ±  2.0  %  as  previously  obtained  from  pCT  analysis  (Figure  5.3  (e)), 
which  is  in  good  agreement  with  that  measured  experimentally. 

5.5  Concluding  Remarks 

This  study  details  (1)  the  mathematical  foundation  for  defining  the  dynamic 
evolution  of  an  optically  suspended  particle  (as  detailed  fully  in  Appendix  A),  and  (2) 
the  experimental  evaluation  of  an  ultrasonic  device  for  the  manufacture  of  short 
fibre  composites.  These  are  discussed  in-turn  below. 

The  mathematical  foundation  illustrated  that  individual  control  can  be 
obtained  for  every  reinforcement  entity,  thus  allowing  highly  complex  fibre 
constructs  to  be  formed.  While  particles  cannot  be  treated  independently  in 
ultrasonic  fields,  these  appear  to  be  more  viable  for  applications  on  the  macroscale. 
A  study  using  water  as  the  host  medium,  indicated  good  agreement  between  the 
mathematical  model  and  the  experimental  studies.  Flowever,  it  should  be  noted  that 
no  exact  prediction  can  be  made  of  the  total  number  of  particles  that  are  drawn 
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towards  a  single  node  but,  to  gain  some  indication,  it  may  be  possible  to  estimate 
the  volume  of  low  pressure  and  compare  this  with  the  particles'  dimensions.  In  the 
case  of  larger  particles,  some  form  of  'bridging'  across  multiple  nodes  is  to  be 
expected,  most  likely  following  the  path  of  least  resistance. 

In  the  second  half  of  the  study,  an  ultrasonic  particle  manipulation  technique 
was  developed  and  employed  to  manufacture  thin,  single  layers  of  short  glass  fibre 
reinforced  polymer  composites.  To  achieve  this,  a  new  type  of  ultrasonic  device  was 
developed,  separating  the  acoustic  system  from  the  resin  cavity  to  allow  for  the  easy 
manufacture  of  multiple  samples.  Furthermore,  an  effective  method  has  been 
outlined  to  manufacture  and  mechanically  characterise  these  thin  laminar  two-phase 
materials.  With  an  8  %  difference  in  stiffness  between  the  two  principal  directions, 
anisotropy  was  demonstrated  for  unidirectionally  reinforced  discontinuous  fibre 
composites  under  uniaxial  tensile  loading.  A  43  %  improvement  in  strength  could  be 
observed  for  samples  tested  along  the  direction  of  fibre  reinforcement,  over  those 
strained  in  the  x-direction,  despite  the  relatively  low  volume  percentage  of  the 
reinforcement  phase. 
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6.  Research  Summary  and  Future  Work 

This  report  details  three  specific  investigations  into  the  generation  of  both 
complex  hierarchical  fibre  architecture  and  controlled  ultrasonic  field  effect 
techniques  (at  different  length  scales)  for  additive  layer  manufacturing  of  biologically 
inspired  short  fibre  reinforced  composites. 

The  work  has  demonstrated  the  successful  formulation  of  Zinc  oxide  (ZnO) 
nanorods  grown  in  aqueous  solutions  onto  the  surface  of  E-glass  short  fibres 
(approximately  3  mm  in  length),  and  a  continuous  borosilicate  glass  fibre  60  mm  in 
length.  The  nanorods  do  not  compromise  the  mechanical  properties  of  the  fibres 
and  furthermore  the  TEOS  coating  could  'heal'  defects  and  cause  a  slight  increase  in 
fibre  cross  sectional  area.  In  the  second  half  of  the  study,  both  theoretical  analysis 
and  experimental  development  were  used  to  evaluate  the  potential  of  field  gradient 
technologies  at  generating  (1)  micron  length  scale  aligned  reinforcement  and  (2) 
nano  length  scale  aligned  composite  materials,  with  regard  to  rapid  prototyping 
manufacture  of  complex  fibre  constructs. 

To  achieve  this,  a  new  type  of  ultrasonic  device  was  developed,  separating  the 
acoustic  system  from  the  resin  cavity  to  allow  for  the  easy  manufacture  of  multiple 
samples.  Furthermore,  an  effective  method  has  been  outlined  to  manufacture  and 
mechanically  characterise  the  micron  length  scale  reinforced  thin  laminar 
composites.  With  an  8  %  difference  in  stiffness  between  the  two  principal  directions, 
anisotropy  was  demonstrated  for  unidirectionally  reinforced  discontinuous  fibre 
composites  under  uniaxial  tensile  loading.  A  43  %  improvement  in  strength  was 
observed  for  samples  tested  parallel  to  the  direction  of  the  fibre  reinforcement  over 
those  strained  normal  to  the  fibre  direction,  despite  the  relatively  low  volume 
percentage  of  the  reinforcement  phase.  This  technique  shows  great  potential  for  the 
low  cost  instantaneous  alignment  of  structural  reinforcement  to  generate  the  light¬ 
weight  high  performance  structures  required  for  future  military  requirements. 

The  orientation  and  distribution  of  reinforcing  entities  in  engineering 
composites  is  key  to  enabling  structural  efficiency,  yet  the  architecture  remains 
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simplistic  when  compared  to  the  distinctive  and  unique  hierarchies  found  in  Nature. 
These  biological  'composite'  materials  achieve  such  configurations  by  accurately 
controlling  the  orientation  of  anisotropic  nano-  and  micro-sized  'building  blocks', 
thereby  reinforcing  the  material  in  specific  directions  to  carry  the  multidirectional 
external  loads  at  different  length  scales.  Capturing  the  design  principles  underlying 
the  exquisite  architecture  of  such  biological  materials  will  overcome  many  of  the 
mechanical  limitations  of  current  engineering  composites. 

This  research  has  investigated  one  possible  route  to  the  realisation  of 
advanced  structural  metamaterials  (i.e.  materials  that  today  do  not  exist  in  either 
engineered  or  biological  forms)  with  controlled  orientation  as  a  'bottom-up' 
manufacturing  technique  for  integration  within  AM.  This  study  has  clearly  illustrated 
the  potential  for  controlled  instantaneous  alignment  of  structural  reinforcement. 
Through  the  application  of  novel  multifunctional  fibres  and  matrix  materials  with 
field  effect  techniques  (ultrasonic/magnetic)  new  structural  designs  will  be  realised 
that  are  hierarchical,  multi-material,  multi-functional,  and  multi-scale  offering 
reduced  weight,  improved  reliability  and  affordability. 

To  achieve  this  'future',  active  advanced  multifunctional  materials,  exhibiting 
programmed  intelligence  in  complex  3D  architectures,  should,  and  could,  be 
developed  via  creative  manufacturing.  This  will  require  (1)  new  hierarchical  fibrous 
metamaterials,  (2)  the  on-going  development  of  field  effect  techniques  and  (3)  the 
development  of  novel  smart  matrices.  One  possible  research  programme  to  develop 
these  materials  of  the  future  is  illustrated  in  Figure  6.1. 
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Figure  6.1:  Overview  of  future  research  concepts:  Left  Column-  WP1:  Top , 
Hierarchical  reinforcement  (ZnO  nano  rods  on  glass  fibres);  Mid,  ultrasonic  field  effect 
prediction;  Lower,  Janus  particles  magnetically  aligned;  Centre  Column-  WP2:  Top, 
acoustic  assembly  pattern  prediction  for  radial  transducers;  Lower,  Ultimaker 
desktop  3D  printer;  Right  Column-  WP3:  Top,  biological  functionally  graded 
architecture;  Lower,  self-folding  architecture. 

These  new  modes  of  assembly,  i.e.  manufacturing  as  a  growth  process,  will 
rely  on  smarter  materials,  not  machines  of  increasing  complexity.  Figure  6.1  shows 
the  three  key  overriding  elements  of  a  future  programme  directed  towards  the 
development  of  metamaterials  for  4D  morphogenesis;  namely  synthetic  processing 
routes  based  on  'complex  hierarchical  metamaterials',  'mechanically  coupled  field- 
effect  assembly'  and  '3D  Additive  Manufacturing'.  It  is  proposed  that  through  the 
use  of  planar  patterning  of  microparticles  and/or  nanoparticles,  it  will  be  possible  to 
achieve  programmable  and  reconfigurable  architectures  through  the  field  effect 
formation  and  control  of  'acoustic  or  magnetic  pixels'.  If  successfully  developed  and 
exploited,  this  would  generate  a  unique  method  for  reducing  lead-times  (through 
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rapid  design  to  production)  and  add  significant  value  across  all  engineering 
disciplines  through  smarter  materials  and  structures,  not  just  smarter  machines. 
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Appendix  A:  Optical  levitation  of  cylindrical  parameters 

A  detailed  derivation  of  the  force  equations,  mathematically  describing  the 
behaviour  of  small  levitated  cylinders  in  a  beam  of  light,  is  presented  below.  Whilst 
closely  following  theoretical  interpretations  by  Gauthier  and  Wallace  [9],  the  present 
approach  embraces  some  added  geometrical  complexity.  The  model  is  considered 
valid  for  all  systems  satisfying  the  condition  d  >  10X0,  where  d  is  the  diameter  of  the 
particle  and  760  the  wavelength  of  electromagnetic  radiation  in  free  space.  A  ray 
diagram  defining  any  geometrical  parameters  is  presented  in  Figure  A.l.  In  addition. 
Figure  A. 2  visually  presents  the  constant  force  fields  that  surround  the  levitated 
particle. 


Figure  A.l:  Schematic  representation  of  axial  force  field  contributions.  FP,  FG,  and  FR  indicate 
the  forces  due  to  photon  scattering,  gravity,  and  drag,  respectively. 
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Appendix  A.l.  Scattering  force 

Mathematically,  the  mechanical  response  of  an  object,  situated  in  a  beam  of  light, 
may  be  ex-  pressed  in  terms  of  Newton's  second  law: 


where  F  is  the  force  acting  on  the  object,  dp  the  object's  overall  momentum  change 
due  to  all  pho-  ton  interactions,  and  dt  the  short  interval  of  time  over  which  dp  is 
measured. 

The  momentum,  p  of  a  single  photon  was  formulated  by  de  Broglie,  who  related 
particle-  and  wave-like  properties  in  the  equation 

h 

p  -  —n 
A0 

where  h  is  Planck's  constant,  77A0  the  wavelength  of  light  in  vacuum,  and  n  the 
refractive  index  associated  with  the  lightwave's  propagation  in  the  medium. 

(A)  The  total  axial  force,  Fz  acting  on  the  particle  may  be  expressed  as  the  sum  of 
four  in-  dependent  contributions, 

Fz  =  Flrz  +  Fltz  +  F2tz  +  F2tz. 

While  Firz  accounts  for  the  reflection  of  photons  at  the  cylinder's  lower  surface,  Fltz 
signifies  the  resultant  force  component  associated  with  refraction;  similarly,  F2rz  and 
F2tz  represent  the  reflection  and  refraction  contributions  at  the  upper  surface.  To 
further  define  the  terms  above,  a  mathematical  formulation  must  first  be  found 
describing  the  momentum  transfer  to  the  particle  upon  re-  flection  and/or  refraction 
of  a  single  photon  at  either  surface.  Using  geometrical  arguments  from  Figure  A. 2, 
the  following  relations  may  be  postulated: 

kpirz  =y-«o  [1  +  cos  (20! )] , 

Ao 

Apifz  =  —  [n0  -  ns  cos  (91  -  02)] , 

Ao 

h 

^P2rz  =  T~ns  [COS  (0!  -  02 ) 

Ao 

+  cos  (302  —  0i)]  , 

fl 

Ap2tz  =  T~  {ns  COS  (01  -  02 ) 

Ao 

-no  cos  [2  (01  -  02)]}  • 
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Figure  A2:  Ray  diagram  schematically  illustrating  the  optical  levitation  of  cylindrical  objects 
inclined  at  an  angle  of  (a)  0  rad,  (b)  rad,  and  (c)  tt/2  rad. 

As  before,  the  indices  z,  1,  2,  r,  and  t  stand  for  axial  direction,  lower  surface,  upper 
surface,  re- flection,  and  refraction,  respectively. 


(B)  The  number  of  photons  acting  on  the  cylinder's  outer  surface  is  directly 
dependent  on  the  physical  characteristics  of  the  light  beam.  Mathematically,  the 
relation  between  these  quan-  tities  is  expressed  in  terms  of  the  beam  intensity 
profile,  I.  For  the  purposes  of  this  study,  the  beam  intensity  profile  may  be 
approximated  as  that  of  the  first  order  Gaussian  mode 

—2  P2' 

W  (z)2 


Kp>z)  = 


2  P 


TtW(z): 


■exp 


where  p  is  the  radial  distance  from  the  beam's  central  axis  and  P  is  the  total  power 
associated  with  the  beam.  \N(z)  describes  the  beam's  width  as  a  function  of  distance, 
z  away  from  the  beam's  minimum  waist  location: 


W(z)  =  W0 


1 

2 


where  z  is  taken  to  be  positive  in  the  direction  of  light  propagation.  W0  is  given  by 


Wo  = 


1 

2 


where  z0  specifies  the  Rayleigh  length  or  the  position  along  the  beam's  axis  at  which 

I  (P/  z)  =  lmax/2. 


(C)  The  power,  P  in  the  light  beam  may  be  written  as 


P  ^av  Ep  h 


with 
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^civ 


N 

dt 


and 


he 

Eph  =  Y 

Here,  N  represents  the  total  number  of  photons,  while  c  stands  for  the  speed  of  light 
in  free  space.  79<t>av  and  Eph  refer  to  the  average  photon  flux  and  the  energy  carried 
by  a  single  photon,  respectively. 

(D)  The  fraction  of  photons  reflected  and  transmitted  at  the  interface  between 
materials  are  given  by  Fresnel's  equations.  The  associated  power  coefficients  may,  in 
the  present  case,  be  written  as: 

I  12  _  9l 

111  W2  +  Q3}2 


with 


Qi  =  ( n0ns)2[cos 2(91)  -  cos2(92)]2 
Qz  =  n0  ns[cos2(G1)  +  cos2(92 )] 
Qs  =  (n02  +  n2)  cos(0i)  cos{92 ) 


and 


\i  |2  1  I  i2 

\ti\  =  1  —  |r*i| 

where  |r|2  and  l^l2  represent  the  total  power  reflection  and  transmission 
coefficients ,  respectively. 

Using  the  results  from  (A)  and  (D),  the  total  change  of  momentum  of  the  particle  in 
the  axial  direction  due  to  its  interaction  with  a  single  photon  can  be  written  in  the 
form: 


dpz  =  \n\2  Apirz  +  \ti\2  Apitz  +  | r  1 1 2 
x  \t-i_l2  Ap2rz  +  |fi|2  |fi|2  Ap2tz 

At  the  same  time,  the  net  axial  force  acting  on  the  particle  may  be  expressed  in 
terms  of  three  fundamental  interactions,  namely  photon  scattering,  gravity,  and 
drag: 


axial  net 


-  F7  -  Fn  -  Fa 
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While  the  subtractive  terms  above  may  readily  be  postulated  as 

Fg  =  v(ps-po)g 


and 


1  _ 

Fg  =  ~p0v2CDA 

the  scattering  component,  Fz  may  now  be  defined  as  the  linear  combination  of  four 
independent  surface  integrals,  which  are  given  by: 


0  L 


F'lrz  ~  I  I  —  [1  +  COS  (2#i)] 

—7 r  0 

x  I  (p,  z )  |ri|2  Rdsdcft, 

o  L 


F\tz  — 


1 

c 


[n0  —  ns  cos  ( 6i  — 


—7 T  0 

x  /  (p,  z)  |fi|2  Rds  dq b, 

o  L 


e2)} 


F2rz  =  f  f  ^  [cos  {e^e2) 

—  71  0 

+  cos  (36>2  —  9i)\  I  ( p,z )  |t i 
x  |fi|2  Rdsdi f>, 


and 


0  L 

F2tz  =  J I  ~ins  cos  _  ^2) 

—7 r  0 

-n0  cos  [2  (6>!  -  62)]}  I  (p,  z) 
x  l^l2  |2  R  ds  dcj). 


R,  ds,  and  d(p80  denote  a  transformation  to  cylindrical  coordinates. 
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Appendix  A. 2.  Gradient  force 

As  with  the  axial  component,  the  total  radial  force,  Fr  acting  on  the  particle  may  be 
expressed  as 

Ft  =  F1tt  +  Fltr  +  F2tt  +  F2tT. 

Using  Figure  A2,  a  subsequent  geometrical  analysis  yields 

Aplrr  =  —  —n0sin  (2 9i)  cos  (ip) , 

Ao 

Apitr  =T -n8  sin  (6>i  -  d2)  cos  (ip) , 

Ao 

Ap2rr  =F-ns  [sin  (30 2  -  6*1 ) 

Ao 

-  sin  (6*1  -  9 2)]  cos  (ip) , 

Ap2tr  =f-  {no  sin  [2  (91  -  6>2)] 

Ao 

— ns  sin  (9 1  —  02)}  cos  (ip) . 

A  new  angle,  MJ  is  introduced  to  map  the  momentum's  radial  direction  onto  the  unit 
vector,  f  defined  in  Figure  7.  Analogous  to  previous  argumentation, 

dPr  =  \ri\2  Aplrr  +  |ti|2  Apltr  +  | r 1 1 2 
X  \ti\2  Ap2rr  +  \tl\2  \tl\2  Ap2tr. 

Not  being  affected  by  gravity,  the  net  radial  force  is  simply  given  by 

Fradial  net  FT  F^ 

Again,  the  scattering  component,  FT  is  obtained  in  terms  of  independent  surface 
integrals  of  polar  form: 
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0  L 


F\rr 


—  sin  (26 1)  cos  (ip) 
c 


—7 r  0 

x  /  (p,  z)  |ri|*  Rds  dcp, 

o  L 


F Hr  =  I  I  ~  Sin  (^1  _  O'l)  COS  W 
—7 r  0 

x  I  (p,  z)  |ti|2  Rds  dtp, 

o  L 


F‘2rr  — 


Ox) 


—  7T  0 

—  sin  (6*1  —  62)]  cos  ( ip ) 
x  I  (p,  z)  | ri  | 2 1 | 2  Rds  d(f>, 


0  L 

F2tr  =  j  j  ~  {n°  Sil1  t2  (01  _  02)] 

—7 r  0 

— sin  (0!  —  02)}  cos  (ip) 
x  /  (p,  z)  |tx|2  |tx | 2  dcj). 


It  should  be  noted  that  p  —  p(9t)  if  the  particle  is  centred  in  the  Gaussian  beam,  but 
that  p  —  p(91,xp)  in  general. 
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Appendix  A. 3.  Torque 

Given  the  particle's  non-spherical  geometry,  a  torque,  T  is  generated  causing  the 
rotation  of  the  cylinder.  To  calculate  this  torque,  the  mathematical  equation 

T  =  Y  X  clF axial  I  Y  X  clF radial 

Js  net  Js  net 

may  be  used.  In  the  above  expression,  5  describes  the  area  over  which  momentum  is 
transferred,  r  defines  a  central  reference  point, 

I(p,z) 

dF  axiai  =  —= - dpzdAez, 

net  h/ph 


and 


dF 


radial  — 
net 


I^ldprdAer. 

Eph 


The  differential  element,  dA  is  representative  of  a  small  surface  element  along  the 
object's  bottom  face.  In  case  the  light  beam's  central  axis  coincides  with  either  of  the 
particle's  axes  of  symmetry,  the  torque  will  vanish. 
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Appendix  B:  Optical  levitation  of  cylindrical:  MATLAB 
development  and  experimental  validation 

Appendix  B.l.  MATLAB  source  code 

MATLAB  source  code  was  developed  to  calculate  the  dynamic  behaviour  of  a 
cylindrically  shaped  object  under  the  influence  of  an  intense  beam  of  light.  Both  the 
changes  in  linear  and  angular  momentum  due  to  light-matter  interactions  between 
the  beam's  photonic  particles  and  the  cylinder's  outer  surface  are  determined  within 
the  source  code.  Furthermore,  any  forces  acting  on  the  levitated  object  due  to 
changes  in  momentum  are  evaluated.  Following  these  computations,  the  progressive 
dynamic  behaviour  of  the  overall  system  is  modelled.  A  number  of  assumptions  are 
made  to  more  readily  enable  coding  of  the  present  simulation,  i.e. 

All  particle  reinforcements  are  rigid  and  cylindrical  in  shape 
The  reinforcement  phase  is  homogeneous 

For  a  given  time  interval,  dt,  the  total  displacement  of  each  fibre  can  be 
expressed  as  the  superposition  of  displacements  due  to  linear  and  angular 
momentum  changes. 

All  photon/body  collisions  are  perfectly  elastic 

Reinforcements  are  free  to  translate  but  may  only  rotate  about  the  cylinder's 
centre  of  mass 

Using  the  optical  field  gradient  developed  by  the  piezoceramic  crystal  array, 
consisting  of  three  sets  of  opposing  ultrasonic  drivers  (T1  -  T6)  (as  illustrated  in 
Figure  3),  the  MATLAB  source  is  able  to  predict  the  resultant  acoustic  field 
depending  upon  the  source  frequency. 
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Figure  Bl.l:  Device  schematic.  T1  -  T6  mark  six  identical  piezoelectric  transducer 
elements.  In  addition,  a  small  piezoceramic  plate  (T7)  provides  particle  levitation 
capabilities.  The  maximum  trapping  area  was  restricted  to  0.225  mm2  . 

Representative  predictions  and  corresponding  experimental  results  (viewed 
as  micrographs)  for  the  interference  of  plane  counter-propagating  travelling  waves 
at  2.300  MHz  for  a  single  set  of  opposing  transducers  pair  (T3  and  T6)  and  two  sets 
of  opposing  transducers  pairs  (Tl,  T4  &  T3,  T6)  are  presented  in  Figure  B1.2.  The 
colour  scale  indicates  absolute  pressure.  The  total  amplitude  is  normalised  to  the 
maximum  attainable  pressure  for  the  present  device  configuration  (see  Figure  Bl.l). 
The  x-  and  y-dimensions  are  measured  in  [pm]. 

Figure  B1.3  (top-left)  indicates  an  alternative  acoustic  pressure  amplitude 
interference  pattern  when  the  T2  &  T5  transducers  emit  a  set  of  plane  counter- 
propagating  travelling  waves  at  2.300  MHz.  The  matching  experimental  results  are 
provided  (top-right).  Figure  B1.3  (bottom)  shows  the  corresponding  result  for  the 
prediction  and  experimental  results  for  T2,  T5  and  T3,  T6  matched  pairs  with  a  plane 
counter-propagating  travelling  waves  at  2.300  MHz 
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Figure  B1.2:  Prediction  and  experimental  results  for  the  interference  of  two  plane 
counter-propagating  travelling  waves  at  2.300MHz.  (Top  Left)  Acoustic  pressure 
amplitude  for  a  single  set  of  opposing  transducers  pair  (T3  and  T6)  with  matching 
experimental  result  (Top  Right).  (Bottom  Left)  Acoustic  pressure  amplitude  for  two 
sets  of  opposing  transducers  pairs  (Tl,  T4  &  T3,  T6)  with  matching  experimental 
result  (Bottom  Right) 
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Figure  B1.3:  Prediction  and  experimental  results  for  the  interference  of  two  plane 
counter-propagating  travelling  waves  at  2.300MHz.  (Top  Left)  Acoustic  pressure 
amplitude  for  a  single  set  of  opposing  transducers  pair  (T2  and  T5)  with  matching 
experimental  result  (Top  Right).  (Bottom  Left)  Acoustic  pressure  amplitude  for  two 
sets  of  opposing  transducers  pairs  (T2,  T5  &  T3,  T6)  with  matching  experimental 
result  (Bottom  Right) 


Both  figures  show  good  agreement  between  the  model  and  the  experimental 
studies.  However,  it  should  be  noted  that  no  exact  prediction  can  be  made  of  the 
total  number  of  particles  that  are  drawn  towards  a  single  node  but,  to  gain  some 
indication,  it  may  be  possible  to  estimate  the  volume  of  low  pressure  and  compare 
this  with  the  particles'  dimensions.  In  the  case  of  larger  particles,  some  form  of 
'bridging'  across  multiple  nodes  is  to  be  expected,  most  likely  following  the  path  of 
least  resistance. 
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